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PREFACE 


In the second decade of the 21st century, NASA plans additional science orbiters, rovers, and 
landers. One proposal is for a Mars sample return mission that would use robotic systems and a 
Mars ascent rocket to collect and send samples of Martian rocks, soils, and atmosphere to Earth for 
detailed chemical and physical analysis. Researchers on Earth could measure chemical and physical 
characteristics much more precisely than they could by remote control. 

Jet Propulsion Laboratory, Mars Exploration Website, 2006 


There has been interest in obtaining gas, dust, rocks, or soil samples from comets or nearby planets 
since the beginning of the space age in the 1950s. Flying a spacecraft to the destination and returning 
a sample in an Earth entry vehicle can accomplish this. The return spacecraft is called a sample return 
vehicle (SRV). Detailed studies of Martian samples can only be conducted on Earth. 

In the past decade there have been many studies investigating the feasibility and practicalities of 
sample return missions to Mars; see Wercinski, 1996; Mitcheltree et ah, 1998; and Desai et ah, 2000. 
Samples collected on Mars will be sealed in a special capsule to prevent contamination. An ascent 
vehicle will launch the capsule from the surface of Mars. There are several ways to return the sample 
to Earth. One way is to use a Mars orbiter to capture the sample container and then begin an 
interplanetary trajectory to return the sample to Earth. 

This document is concerned with the Earth entry environment and the Earth entry portion of a sample 
return mission to the planet Mars. The document is intended to introduce the subject of atmospheric 
entry to engineers and scientists who do not have strong backgrounds in aerodynamics, 
aero thermodynamics and flight mechanics. Much of this material can be applied to entry into other 
planetary atmospheres by using the appropriate atmospheric and gravity models. This document 
discusses topics in atmospheric entry and presents examples in related disciplines. References are 
cited for further infonnation. The document is not intended to be comprehensive and some important 
topics are omitted (e.g., chemically reacting gases and nonequilibrium flows). There is an extremely 
large number of engineering and scientific books, journal articles, conference papers, and 
government, industry and university reports on the subject of atmospheric entry (see National 
Aeronautics and Space Administration (NASA), American Institute of Aeronautics and Astronautics 
(AIAA), Department of Defense, and commercial online databases for several thousand citations on 
atmospheric entry and/or re-entry vehicles). 

The topics considered in this document include basic principles of physics (fluid mechanics, rigid 
body dynamics, and heat transfer), chemistry, and engineering. The following subjects are discussed: 

Fluid mechanics - Aerodynamics, compressible fluids, shock waves, continuum flow 
(incompressible to hypersonic), hypervelocity flow (high temperature gas dynamics), boundary 
layers and viscosity, and free molecular flow. 

Earth ’s atmosphere and gravity - Earth model atmosphere and acceleration due to gravity. 
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Heat and mass transfer - Forced convection, radiation, and ablation. 

Dynamics - Rigid body dynamics, and static and dynamic stability. 

Flight and ground test simulation - Flight simulation and ground test facility capabilities. 

Numerical simulation - A brief description of computational fluid dynamics (CFD) and trajectory 
simulation tools. 

For definitions of unfamiliar terms see Definition of Technical Terms for Aerospace Use (Allen, 
1965), or the current online version (Glover, 2006). 

A brief historical note may help to introduce the material that follows. In October 1957 the Soviet 
Union launched the first artificial satellite, Sputnik I, and the space age began. In 1958 the U.S. 
Congress passed the Space Act creating the National Aeronautics and Space Administration (NASA). 
In 1961, President John Kennedy announced that the U.S. would land a man on the moon and return 
him safely to Earth before the end of the decade, a feat accomplished by the Apollo 1 1 mission in 
July 1969. 

In the 1960s the U.S. and the Soviet Union developed intercontinental ballistic missiles (ICBMs) 
capable of flying more than 10,000 km in 30 minutes. The U.S. solid propellant Minuteman II was 
operational in the 1960s and was capable of a range of 1 1,000 km. The Minuteman II weighed 
31,750 kg at liftoff and traveled at 6.7 km/sec at third-stage burnout. The two-stage, liquid propellant 
Titan II, also operational in the 1960s, was capable of a range of 15,000 km and reached an altitude 
of 1500 km. The burnout velocity for the Titan II was greater than 6.7 km/sec and liftoff weight was 
149,700 kg (Blake, 1988). Atlas, Titan and Minuteman ICBMs have been used as launch vehicles for 
satellites and spacecraft. By the 1970s Soviet ICBMs were capable of 16,000-km ranges with liquid 
propellant rockets. 

Prior to 1958, the National Advisory Committee for Aeronautics (NACA), the predecessor of NASA, 
was involved in atmospheric-entry research. NACA engineers and scientists experimentally and 
analytically studied supersonic and hypersonic flows. The “blunt body concept” to reduce 
aerodynamic entry heating was formulated by NACA’s H. J. Allen in the 1950s (Allen and Eggers, 
1958). Allen demonstrated that blunt, high-drag bodies are superior to pointed, slender bodies in 
atmospheric entry. This is due to the fact that the stagnation-point, convective heat flux decreases as 
bluntness increases. The higher drag of the blunt body also tends to reduce the convective heat flux 
by causing the body to decelerate at higher altitudes (see sec. 4.2). The Mercury, Gemini and Apollo 
entry capsules were applications of the blunt body concept. The early days of the space age and the 
formation of NASA are described in a recent NASA monograph (Portree, 1998). 

Lander spacecraft are designed to reach the surface of a planet and survive long enough to record and 
transmit scientific data. The first successful planetary landing was by the Soviet Union on Venus in 
1966. The Soviet Venera landers survived the harsh conditions on Venus while carrying out chemical 
composition analyses of the rocks and relaying color images. The Soviets also landed spacecraft on 
Mars in 1974. The United States' Mars Viking landers in 1976 and the earlier Surveyor landers on the 
moon carried out similar experiments. NASA participated extensively in the design and development 
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of many successful planetary probes including the Viking Mars lander, Pioneer Venus, Galileo 
(Jupiter arrival in 1995) and Mars Pathfinder in 1997. 

The Space Transportation System (STS), or Space Shuttle Orbiter, was the subject of extensive 
ground testing and numerical simulation. An extremely important contribution to the Orbiter’ s 
success was the NASA Ames-Lockheed development of a revolutionary thermal protection system 
(TPS). Lockheed furnished lightweight ceramic tiles and Ames provided tile coatings and flexible 
blankets to shield the vehicle and crew from the intense entry heat loads. Similar TPS technology 
developments made possible many successful interplanetary missions. There have also been 
remarkable advances in computational fluid mechanics and in information technology in the past 
decade. These facts are mentioned because they will have a significant influence on future space 
exploration. 

This document is an attempt to describe some of the basic concepts in aerodynamics and 

aero thermodynamics that relate to the entry of a Mars sample return vehicle (SRV) into the Earth's 

atmosphere. 

A personal note - When the first Sputnik was launched I was a graduate student in engineering. In 
the years that followed there were many changes in the engineering curriculum. Engineering colleges 
began to offer courses in celestial mechanics and “space technology.” University departments of 
aeronautics and mechanical engineering added aerospace engineering and astronautics to the 
curriculum. More emphasis was placed on hypervelocity flight, rarefied (low density or free 
molecular) flow and high temperature gas dynamics. Journals, conference papers, courses and 
textbooks were devoted to these subjects. Universities and industry began to place greater emphasis 
on collaboration between scientists and engineers. 

Note on Internet resources - There is a great deal of useful information on aeronautics and space on 
the Internet (World Wide Web). Unfortunately, websites change frequently and often disappear. 
Several websites are mentioned in this document and listed in the References and Bibliography 
sections. Most of the sites listed are those of professional societies and other well-established 
organizations. 
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NASA Photo No. A-22664 (1957) 

II. Julian Allen is best known for the "Blunt Body Theory, ” an aerodynamic design technique for 
alleviating the severe re-entry heating, which was delaying the development of ballistic missiles. His 
findings revolutionized the fundamental design of ballistic missile re-entry vehicles. Subsequent 
research led to the use of blunt shapes for ballistic missiles and spacecraft re-entering the Earth's 
atmosphere. Blunt, ablative heat shields were used to protect the Mercury, Gemini and Apollo 
astronauts as their space capsules returned to Earth. 


xii 


1. INTRODUCTION 


1.1 Scope 

This document is intended to provide an introduction to Earth atmospheric entry for the 
nonspecialist. However, a basic knowledge of engineering or physics is required. The purpose of this 
work is to describe various aerodynamic, aerothermodynamic, and flight-dynamic aspects of the 
high-speed entry of sample return vehicles (SRVs) into the Earth’s atmosphere. SRVs returning from 
missions to Mars would have entry velocities of 1 1 to 15 km/s. 1 Flight velocities of this magnitude 
are often referred to as hypervelocity rather than hypersonic. In hypervelocity flight high- temperature 
chemical reactions significantly affect the airflow around an entry vehicle. At altitudes above 80 km 
the density of the atmosphere is so low that the average distance an air molecule travels between 
successive collisions with other air molecules is large compared to the size of the SRV. This 
condition is referred to as the free molecular flow regime in contrast to the continuum flow regime at 
lower altitudes. 

Aerodynamic and gravitational forces strongly influence SRV entry dynamics and heat transfer 
characteristics. Knowledge of atmospheric properties is necessary to detennine aerodynamic and 
heating loads. Aerodynamic forces and heat transfer characteristics are functions of atmospheric 
density and vehicle speed. A model Earth atmosphere is discussed briefly in section 2. 

The design an SRV for a specific mission requires knowledge of the characteristics of the flow 
around the vehicle. Detailed flow field information may not be necessary if experimental data are 
available for similar vehicles and flight conditions. Computational fluid dynamics (CFD) 2 solutions 
that have been validated can also be used in the vehicle design process. 

Preliminary trajectory calculations can be used to predict the flow field around a proposed vehicle. 
Knowledge of the vehicle flow field at several points along the trajectory enables computation of 
vehicle surface pressures, aerodynamic forces, and heat transfer characteristics. Mach number and 
Reynolds number 3 histories are also necessary to: a) characterize the vehicle boundary layer (laminar, 
transitional or turbulent); b) provide initial conditions for CFD solutions; and c) establish appropriate 
requirements for ground testing. 

In the preliminary design phases the mission requirements establish a vehicle shape based on payload 
volume and aerothermodynamic limitations for a proposed trajectory. Costs, schedules, launch 
vehicle availability, reliability and recovery considerations are factored into the design. Development 
testing and analysis are performed as necessary. As the design progresses more detailed trajectory 
information becomes available to ensure that stability, guidance, communications, navigation, and 


1 At an altitude of 150 km (93.2 mi) a spacecraft or satellite in a circular Earth orbit travels at a speed of 7.8 km/s 
(25,600 ft/s). In order for a spacecraft to escape the Earth’s gravitational field the spacecraft must reach a speed of 
approximately 1 1.2 km/s; this is referred to as the Earth escape velocity. 

2 Computational fluid dynamics is the discipline that uses high-speed digital computers to solve the equations of flow 
around or within various bodies. See section 1.7 for a brief description of computational fluid dynamics. 

3 Mach number is the ratio of the vehicle speed to the speed of sound in the surrounding fluid. It is a measure of the 
compressibility of the fluid. Reynolds number is the ratio of inertia (or momentum) forces to viscous forces in a fluid 
flow. It is a function of the speed and size of a body and the viscosity of the fluid. Mach number and Reynolds number 
are defined in section 1.3.1. 
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control requirements are satisfied. Validation and qualification tests are performed and the design is 
baselined. Note that aerothermodynamic design considerations and trajectory calculations are critical 
in the early and later design phases. Other aspects of spacecraft systems design are discussed in 
section 3. 


1.2 Basic Considerations 

Sample return vehicles are designed to ensure that the vehicle payload will survive atmospheric 
entry. The following conditions must be carefully assessed in the design of an SRV: a) maximum 
aerodynamic load; b) maximum deceleration load; c) maximum heat flux; and d) total heat load. The 
maximum aerodynamic load usually occurs when dynamic pressure is a maximum (aerodynamic 
forces are directly proportional to dynamic pressure). If the vehicle angle of attack 4 is varying, it is 
necessary to examine the product of dynamic pressure and angle of attack (lift force usually increases 
as angle of attack increases). The maximum aerodynamic load usually occurs when the product of 
dynamic pressure and angle of attack is a maximum. 

Excessive deceleration loads can be decreased by increasing the vehicle lift-drag ratio L/D. Lifting 
bodies can be maneuvered to control the SRV impact location. Controlling lift and bank angle 5 can 
vary both down-range and cross-range impact locations. 

Providing an adequate thermal protection system (TPS) is a major design consideration. The 
structural and heat shield materials selected must be able to withstand the maximum heat flux and the 
total heat load. Greater heat fluxes are encountered in steep entry than for shallow entry angles. 6 
However, the time of flight is usually greater for shallow entry angles, thus increasing the total heat 
load the vehicle must absorb (since total heat load is equal to the integral of heat flux). 7 

At SRV entry speeds greater than 1 1 km/s temperatures greater than 1 1,000 K are encountered in the 
shock layer. (The region between the bow shock and the vehicle surface is called the shock layer; see 
fig. 1 .) At these temperatures, the heat flux to the vehicle due to radiation from shock-layer gases can 
be significant. At lower velocities and temperatures the heat flux is primarily convective. Ceramic 
tiles used on the Space Transportation System (STS) Orbiter and ablating materials 8 used on the 
Apollo command module and planetary entry vehicles protected these vehicles from high 
temperatures and high heat loads. Figure 2 shows shock layer temperature as a function of altitude 
and vehicle velocity for equilibrium air (equilibrium and nonequilibrium flows are discussed in 
sec. 1.3.2). Note that thermodynamic and transport properties (e.g., viscosity and thennal 


4 The angle of attack is the angle between the vehicle axis and the vehicle velocity vector (see fig. 12). 

5 The bank, or roll angle is the angle between a vertical (or horizontal) reference plane and the lift vector. 

6 The flightpath angle is the angle between the horizontal and the vehicle velocity vector; the positive direction is up and 
the flightpath angle at entry is always negative. (A flightpath angle of-20 deg is shallower than a flightpath angle of 
-40 deg.) 

7 In this document, the heat flux q is defined as the heat transfer rate per unit area. The symbol q (or q-dot) is often used 
for heat flux. The heat load Q is the time integral of heat flux. In forced convection heat transfer, the heat flux is 
proportional to the temperature gradient normal to the flow direction at the body surface. The constant of 
proportionality k is the transport property known as the thermal conductivity of the fluid. 

8 Ablating materials are materials used on the surface of an entry vehicles to absorb heat by removal of mass, thus 
blocking the transfer of heat to the vehicle. The total heat load determines the amount of material required to protect the 
vehicle. 
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conductivity) of air are not well understood at temperatures above 9000 K. Other features of 
hypersonic flow around a blunt body are described in the sections that follow. 


1.3 Flow Regimes 

Atmospheric entry trajectories are shown as functions of altitude and velocity in figure 3 for Space 
Shuttle Orbiter, Apollo 4, and for a typical SRV (solid lines). Figure 3 also shows approximate 
boundaries for vibrational excitation, oxygen and nitrogen dissociation, ionization, and 
thermochemical nonequilibrium (dashed lines). At high altitudes free molecular flow effects are 
significant (above 150 km for the Orbiter) and continuum flow (see sec. 1.3.1) approximations are 
not valid. At low altitudes continuum subsonic, transonic, supersonic and hypersonic flows occur. At 
vehicle speeds below about 100 m/s the flow can be considered incompressible (air density is 
assumed constant); water is an example of an incompressible fluid. For each flow regime it is 
necessary to consider the flow field around the vehicle, the vehicle surface pressure distribution, and 
the heat transfer to or from the vehicle. 

1.3.1 Continuum Flow 

For air at standard sea level conditions, the distance a molecule travels between collisions is 
extremely small and the atmosphere appears to be a continuous medium to an object moving through 
it. 9 Continuum flow occurs when the molecular mean free path X is much smaller than a characteristic 
vehicle dimension L; i.e., X « L. Knudsen number Kn is defined as the ratio of A, to L. Knudsen 
number can also be expressed as a function of Mach number and Reynolds number. Continuum flow 
occurs when Kn « 1 and free molecular flow occurs when Kn » 1 . In terms of Mach number M 
and Reynolds number Re, continuum flow occurs when M/(Re) 1/2 < 0.01 and free molecular flow 
occurs when M/Re > 3 (Regan and Anandakrishnan, 1993); see table 1. Free molecular and 
continuum flow regime boundaries are shown as functions of Mach number and Reynolds number in 
figure 4. Mach number is defined below. Reynolds number is defined in section 1.5. 

TABLE 1. FLOW REGIME BOUNDARIES 


Flow Regime 

Anderson (1989) a 

Regan and Anandakrishnan (1993) 

Free molecular 

Kn » 1 

Re < M/3 

Near free molecular 

Kn > 1.0 

— 

Transitional 

1.0 >Kn> 0.03 

M/3 <Re< 10,000 M 2 

Continuum 

Kn < 0.2 

Re > 10,000 M 2 


a Note the overlap of the continuum and transitional flow regimes. Anderson (1989) notes that in this region 
(0.03 < Kn < 0.2) temperature and velocity slip effects (discontinuities) are present at the body surface. 


9 At sea level the molecular mean free path for air is 6.633 x 10 s m; at 100 km altitude the mean free path is 0.141 m (see 
National Oceanic and Atmospheric Administration, U.S. Standard Atmosphere, 1976). 
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A perfect gas is both thermally perfect and calorically perfect. A thermally perfect gas obeys an ideal 
gas equation of state. A calorically perfect gas has constant values of specific heat, independent of 
temperature; i.e., specific heat at constant pressure c p and specific heat at constant volume c v are 
constant. A perfect gas is a gas that: satisfies an equation of state; has specific heats independent of 
temperature; and has internal energy that is a function of temperature alone. A perfect gas is also 
called an ideal gas. 10 (see Allen, 1965 and Glover, 2006). A nonperfect gas is often referred to as a 
real gas. 

The equation of state for a thermally perfect gas is 

p = p R T 


where 


p = pressure 

p = mass density 

T = absolute temperature 

R = gas constant for a specific gas 

R = R / M = 287 J/kg-K for air 

R = universal gas constant = 83 14.34 J/(kmol-K) 

M = molecular weight = 28.9644 kg/kmol (for air at altitudes <85 km). 

Note that the perfect gas equation of state given above is not valid for hypersonic flow (Mach 
numbers above five or six) where the isentropic exponent (or ratio of specific heats) y for air is no 
longer constant. Figure 5 shows the variation of y with shock layer temperature and pressure for air 
(Gordon and McBride, 1994). Only at temperatures below 800 K (or velocities below about 800 m/s) 
can air be assumed to be calorically perfect. 

In continuum flow Mach number M is a very important dimensionless parameter. Mach number is 
defined as the ratio of the vehicle (or freestream) speed V to the local speed of sound c: 

M = V/c 


where 

c = (y R T) ‘ for a thermally perfect gas 
y = isentropic exponent = Cp/c v 
c p = specific heat at constant pressure 
c v = specific heat at constant volume. 


10 The terms perfect fluid and ideal fluid usually refer to an inviscid fluid that may also be incompressible. 


4 



An approximation for the speed of sound in thermally perfect air with y = 1.4 and T < 1000 K is 

c * 20.0 (T) 1/2 m/s . 

Figure 2 shows the effect of freestream, or vehicle, velocity on shock layer temperature. Figure 5 
shows the effect of shock layer temperature and pressure on the isentropic exponent y for air. These 
variations are important considerations in hypersonic flow. 

Mach number is a measure of the effect of the compressibility of a fluid. When Mach number is near 
or equal to zero the fluid is considered incompressible and the density of the fluid is assumed 
constant. Flows with Mach number less than one and without shock waves are classified as subsonic. 
As the Mach number increases, compressibility effects become significant and density can no longer 
be assumed constant (compressible subsonic flow). 

Analytical solutions for incompressible flow around various body shapes can be obtained by solving 
the Laplace equation for the velocity potential (Karamcheti, 1980). Numerical solutions for subsonic 
flow around arbitrary two- and three-dimensional bodies can be obtained using panel methods (Katz 
andPlotkin, 1991). 

In incompressible flow static pressure p can be obtained from Bernoulli’s equation: 

p + 1/2 p V 2 = const 


where 


p = const for M « 1 . 


For compressible flow (p ^ const) 


y/(y - 1) p/p + 1/2 V 2 = const 


It is often useful to determine a pressure coefficient distribution to characterize a flow field. The 
pressure coefficient is defined as 


Cp — (p Poo) / qoo 


where q ,. is the freestream dynamic pressure 

2 

Qoo = 1/2 Poo Voo 


or for a thermally perfect gas 


Qoo 1/2 y Poo Moc . 
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Note that the subscript oo indicates conditions far upstream of the body. The subscript is dropped in 
later sections when it is not necessary to distinguish between freestream and local conditions. 

Similarity rules - Similarity rules can be used to determine lift, drag and pressure coefficients from 
ground test results or from CFD solutions for similar body shapes. That is, the results for one body 
shape can be scaled to another slightly thicker or thinner similar body at a slightly lower or higher 
Mach number. 

For compressible subsonic flow (Mach numbers less than about 0.7), pressure and force coefficients 
on slender bodies at small angles of attack a can be scaled using the Prandtl-Glauert and Goethert 
similarity rules (Liepmann and Roshko, 1957). 

Body 1 has a fineness ratio (body thickness/body length) t = Ti, a freestream Mach number M = Mi, 
and a pressure coefficient C p = C p i; and similarly for Body 2. Assume T| and x 2 « 1 and ai and 
a .2 « 1 radian. 

a) If x 2 / (1 - M 2 2 ) 1/2 = n / (1 - Mi 2 ) 172 , then 


If Mi « 0, then 


C p2 = C p i 


C p2 = C p i = C 


p INCOMP 


b) If x 2 = Ti, then 


C p2 = C p i (1 - Mi 2 ) 172 / (1 - M 2 2 ) 172 


If Mi - 0, then C p i - C p incomp and 


C p2 


C p inCOMp/ (1 - M 2 -) 17 


c) If M 2 = Mi, then 


C p2 C p i t 2 /ti 

The Goethert rule applies to slender, axially-symmetric, as well as two-dimensional, bodies: 

Ifx 2 = Ti(l — Mi 2 ) 172 /(1 -M 2 2 ) 172 then 

C p2 = C p i (1 - Mi 2 ) 172 / (1 - M 2 2 ) 172 

There are also transonic, supersonic and hypersonic similarity rules (Van Dyke, 1951 and Spreiter, 
1982). The supersonic similarity rule is similar to the subsonic rule with (1 - M~) replaced by 
(M 2 - l) 172 . 
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Compressibility and flow regimes - Compressibility is the fluid property that relates changes in 
pressure and specific volume (or fluid density). Liquids are highly incompressible and gasses are 
highly compressible. However, at very low speeds gasses behave much like an incompressible fluid. 
As pressure increases specific volume decreases (and density, the inverse of specific volume, 
increases). 

An incompressible flow is one in which the fluid density remains constant. At low subsonic speeds 
(M « 1) fluids can be considered incompressible. As vehicle speed increases, changes in the density 
of the fluid can no longer be neglected and the flow is referred to as compressible (Glover, 2006). 

As the speed of a vehicle approaches the speed of sound (M = 1 or “Mach one”) shock waves begin 
to fonn and the flow becomes transonic. This is the region of the so-called “sound barrier,” where a 
vehicle experiences a large increase in drag and often experiences unsteady buffeting. 

Compressibility effects become even more significant at Mach numbers above one. Shock waves are 
present and distinct and the flow is called supersonic or hypersonic. In supersonic continuum flow 
the thickness of the shock wave is negligible and flow properties are treated as discontinuities across 
the shock. Shock waves can be attached (pointed body), detached (blunt body), normal, or oblique 
(see sec. 1.5). A pointed body can have a detached shock at low supersonic Mach numbers or at any 
Mach number if the nose semivertex angle is large. A normal shock wave occurs when a supersonic 
flow is brought to rest at a stagnation point 11 on a body surface normal to the flow direction (or a 
stagnation line on a wing-like, two-dimensional body). A normal shock wave can also occur in a 
shock tube when a diaphragm separating high and low pressure gases is removed or ruptured (i.e., as 
the high pressure gas expands into the low pressure region). Attached oblique shock waves occur 
when a supersonic flow impinges on a wedge or an inclined plane (if the flow deflection angle is not 
too large or the Mach number is not too low). 

At Mach numbers above around five, the flow is usually considered hypersonic. (The number five is 
somewhat arbitrary and may depend on the “bluntness” of the body.) Hypersonic flows are 
characterized by shock waves very close to the body surface (thin shock layer), thick boundary 
layers, high temperatures, and aerodynamic coefficients that can be nonlinear functions of angle of 
attack. The “thick” and “thin” characterizations are with respect to a body dimension such as 
diameter or length. Some features of hypersonic flow over a blunt body are shown in figure 1. At the 
rearward facing base of the body the flow separates and creates a region of recirculating flow 
bounded by dividing streamlines. Note that viscous effects are confined to the boundary layer and 
wake regions. 

1.3.2 Hypervelocity Flow 

At temperatures above about 800 K at an air pressure of 1 atm, vibrational excitation of (A and N 2 
molecules becomes significant and specific heats c p and c v are no longer constant. The perfect gas 
equation of state given above continues to apply until dissociation of 0 2 begins at about 2500 K. In 
this temperature range the gas is thermally perfect, but not calorically perfect (Anderson, 1989). At 
temperatures above about 2500 K the gas is a mixture of chemically reacting gases and each 


1 1 A stagnation point is a point in the flow about a body where the fluid particles have zero velocity with respect to the 
body. 
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chemical species obeys the perfect gas equation. For chemically reacting gases it is important to track 
the concentrations of each species. The velocities producing very high temperatures are often referred 
to as hypervelocities rather than hypersonic velocities. 

Dissociation - The dashed lines on the left side of figure 3 indicate that oxygen and nitrogen 
molecules in the air behind (or downstream of) a normal shock wave begin to dissociate to form 
atoms. This is the first significant departure from a calorically perfect gas. For air at a pressure of 
1 atm, O 2 dissociation occurs between 2500 K and 4000 K; N 2 dissociation occurs between 4000 K 
and 9000 K (Anderson, 1989). 

Ionization - To the right of the dashed line on the right side of figure 3 ionization effects become 
important. Electrons leave the oxygen and nitrogen atoms to form ions. This is the cause of the well- 
known communications “blackout” that occurs during re-entry. For air at a pressure of 1 atm, 
ionization of O and N molecules occurs at temperatures above 9000 K (Anderson, 1989). Note that 
NO molecules can become ionized at lower temperatures and pressures. 

Nonequilibrium flow - An equilibrium gas flow is one in which the composition of the gas at any 
point is independent of time and energy is constant along streamlines. A nonequilibrium flow is one 
that is not in thermal equilibrium or chemical equilibrium; this occurs at very high temperatures 
and/or in very high-speed flows (large kinetic energies). 

Thermal, or thermodynamic, equilibrium means that physical properties (e.g., pressure and 
temperature) have ceased to change with time at a given point; or that they are varying so slowly that 
for any point, one can assume thermodynamic equilibrium in some neighborhood about that point. 
Chemical equilibrium means all chemical reactions are in balance and the system does not undergo 
any further change in chemical composition or in the concentrations of the various substances. 

Frozen flow is a special case of nonequilibrium flow. In frozen flow all chemical reactions are 
assumed to have stopped — frozen in time. Air in a frozen state can have properties significantly 
different than equilibrium air at the same temperature and pressure. 

The frequency of molecular collisions decreases as altitude increases; and at altitudes above about 
50 km, the rate of physical and chemical processes must be considered. A finite time is required for 
chemical reactions to reach an equilibrium composition. For a brief interval (a few milliseconds) 
after the temperature increase, the process will be in chemical nonequilibrium. In nonequilibrium 
flows, specie concentrations and thermodynamic properties are functions of time. Equilibrium 
properties for chemically reacting gases can be derived from statistical thermodynamics. 
Computational modeling of nonequilibrium thermodynamics is very difficult and equilibrium 
thennodynamics is often used as a simple approximation to the actual flow. 

Other hyp ervelo city flow effects - Other effects that must be considered in hypervelocity flow 
include: a) the need to dissipate enormous quantities of energy (a very large percentage of the vehicle 
kinetic energy is transferred to the gas in the shock layer); b) very high shock layer temperatures (as 
high as 14,000 K for sample return missions) are sources of heat radiation that can have a significant 
effect on vehicle surface temperatures; and c) foreign gases (from vaporization), solid particles (from 
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spallation 12 due to thennal stresses), and products of ablation produce additional complications in the 
flow field. 

1.3.3 Free Molecular Flow 

At high altitudes collisions between molecules are less frequent and the significant flow parameter is 
the Knudsen number Kn. The Knudsen number is the ratio of the molecular mean free path X and a 
characteristic vehicle dimension L (usually nose radius or body length): 


Kn = X/L 


In general, free molecular flow occurs when Knudsen number is much greater than 1 and continuum 
flow occurs when Knudsen number is much less than 1 (see table 1). In free molecular flow a distinct 
bow shock does not exist in front of a blunt body (Rohsenow et al., 1985). 

In free molecular flow and at altitudes just below the free molecular limit, molecules near the vehicle 
surface may have a mean velocity different from that of the surface. This is referred to as velocity 
slip. If the fluid molecules do not acquire the energy corresponding to the energy of the vehicle 
surface after one collision, there is lack of momentum accommodation. If a temperature discontinuity 
occurs at the vehicle surface, there is a lack of thennal accommodation. 

In free molecular and near free molecular flows, thermal and momentum accommodation 
coefficients 13 are can be determined experimentally and used to predict temperature jump and 
velocity slip at the surface. Thermal accommodation coefficients are functions of gas temperature T, 
surface temperature T w , gas composition, and surface material. For air with T « T w ~ 300 K, the 
thennal accommodation coefficient varies from 0.89 to 0.97 for polished metal surfaces 
(Rohsenow et al., 1985). 

The convective heat flux in free molecular flow can be determined from the kinetic theory of gases or 
obtained from an experimental conelation of Stanton number St and thermal accommodation 
coefficient. The Stanton number is the ratio of the convective heat transfer at the surface to the heat 
transfened by the heat capacity of the gas. 

St = (convective heat transfer coefficient) / p c p V 


where 


convective heat transfer coefficient = q / (Tw - Tp) 


and 


T F = characteristic fluid temperature (usually mean or freestream). 


12 High thermal stresses at the vehicle surface may cause solid particles to spall and enter the flow field. 

13 The accommodation coefficient is the ratio of the average energy actually transferred between a surface and impinging 
gas molecules which are scattered by the surface to the average energy which would theoretically be transferred if the 
impinging molecules reached complete thermal equilibrium with the surface before leaving the surface (Glover, 2006). 
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Free molecular drag coefficients can be detennined from momentum accommodation coefficients. 
Algorithms for calculating free molecular lift and drag coefficients are given in Regan and 
Anandakrishnan, 1993. Direct simulation Monte Carlo 14 (DSMC) computational methods are used to 
determine free molecular and near free molecular flow-field characteristics. 

1.3.4 Transition from Free Molecular to Continuum Flow 

A vehicle leaves the free molecular flow regime as it descends into the sensible atmosphere. Since 
Knudsen number is inversely proportional to a vehicle characteristic length, smaller bodies 
experience free molecular flow at lower altitudes (see table 2). The Space Shuttle Orbiter experiences 
free molecular flow at altitudes above 170 km and continuum flow at altitudes below around 90 km. 


TABLE 2. APPROXIMATE ALTITUDE LIMITS FOR FREE MOLECULAR, TRANSITIONAL, 
AND CONTINUUM FLOW REGIMES, 1976 U.S. STANDARD ATMOSPHERE* 


Characteristic 
body length, 
L, m 

Free molecular flow 
regime, altitude, km 
(Kn> 100) 

Transitional flow regime, 3 
altitude, km 
(1 >Kn>0.03) 

Continuum flow 
regime, altitude, km 
(Kn < 0.01) 

1 

>170 

110-90 

<85 

0.5 

>160 

105-85 

<80 

0.1 

>130 

95-75 

<70 


* National Oceanic and Atmospheric Administration, 1976. 
a Anderson, 1989. 


At altitudes just below the free molecular limit velocity slip and lack of momentum accommodation 
may occur near the vehicle surface. Near the free molecular flow limit there may also be a jump in 
temperature at the vehicle surface. 

At altitudes just below the free molecular flow limit there is a transitional flow regime. (This is not 
the same as the transition from laminar flow to turbulent flow in the boundary layer that occurs as 
Reynolds number increases.) Accurate prediction of aerodynamic coefficients and heat fluxes is 
difficult near the free molecular flow limit and in the transitional flow regime. 

At altitudes just below the transitional flow regime there is a viscous merged layer flow regime. In 
the transitional flow regime shock waves can have appreciable thickness. (In continuum flow the 
thickness of a shock wave is negligible and flow properties are treated as discontinuities at the shock 
location). In the viscous merged layer regime there may be strong interactions between the shock 
wave and the boundary layer. (Regan, 1984 and Anderson, 1989). Approximate flow regime 
boundaries are given in tables 1 and 2 and shown in figure 4. 


14 The direct simulation Monte Carlo (or DSMC) method is used to simulate gas flows in the free molecular flow regime. 
It is a computational technique that models the statistical behavior of the gas molecules. 
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1.4 Shock Waves 


Shock waves are a prominent feature of re-entry flows. In continuum flow the thickness of a shock 
wave is neglected and flow properties are discontinuous at the shock location. A normal shock wave 
is shown schematically in figure 6a. (A nonnal shock wave is perpendicular to the upstream velocity 
vector.) Static pressure p, static density p, and static temperature T increase across the shock wave 
while Mach number M and flow speed V decrease (p 2 > pi, p 2 > pi, T 2 > Ti, M 2 < Mi, V 2 < Vi) 
where the subscripts 1 and 2 indicate conditions just upstream and just downstream of the shock, 
respectively. For a normal shock wave, the downstream Mach number M 2 is always subsonic 
(M 2 < 1). A normal shock wave occurs ahead of a blunt body in supersonic or hypersonic flow and is 
referred to as a detached or bow shock (see fig. 1). The distance between a detached shock and a 
body is called the shock-wave standoff distance. 

An oblique shock wave occurs when the flow direction changes across the shock and the shock wave 
is no longer normal to the freestream direction as shown figure 6b. The component of the velocity 
vector that is normal to an oblique shock wave must be supersonic upstream of the shock and 
subsonic downstream of the shock. Note that the Mach number downstream of an oblique shock 
wave may be supersonic; it is only the velocity component nonnal to the oblique shock wave that is 
subsonic. 15 


1.5 Boundary Layers and Viscous Effects 

Another very important variable in fluid dynamics is the Reynolds number, which is a measure of the 
effect of fluid viscosity. The nondimensional Reynolds number Re is the ratio of inertia forces to 
viscous forces; viscous effects become less important as Reynolds number increases. Reynolds 
number is defined as 


Re = p V L / p 

where 

p = mass density of fluid 
V = speed of flow (or vehicle speed) 

L = characteristic vehicle dimension 
p = absolute viscosity of fluid = p(T) . 

The Sutherland approximation for the absolute viscosity of air is given in section 2. 


15 Equations, tables and charts for flow properties across normal and oblique shock waves are given in NACA Report 
1135 (Ames Aeronautical Laboratory, 1953). Scanned images of NACA reports are available on the NASA Technical 
Reports Server (see References). A website is available to calculate conditions across normal and oblique shock waves 
and other compressible flow properties (Delft Technical University, 2006). 
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Reynolds number can also be expressed in terms of kinematic viscosity v 


Re = V L / v 


where 

v = |u/p. 

In inviscid (nonviscous) flows, viscosity effects are neglected and the fluid velocity at the surface of 
a body can have a finite, nonzero value. In reality, fluid viscosity requires the velocity to be equal to 
zero at the body surface (no-slip condition). At moderately high Reynolds numbers (about 100,000) 
the effect of viscosity is appreciable only in the vicinity of the body. This discovery by the German 
scientist Ludwig Prandtl was one of the most significant in aeronautics and is the basis for boundary 
layer theory. The development of a boundary layer in the flow over a flat plate is shown 
schematically in figure 7. 

In figure 7 the flow along the plate starts out laminar at the leading-edge stagnation point where the 
flow velocity is equal to zero. As the local Reynolds number Re x based on the distance x from the 
leading edge increases, the flow ceases to be laminar and passes through a “transition zone” before 
becoming fully turbulent. 16 Boundary layer transition is affected by several factors including the 
roughness of the body surface and the turbulence in the freestream. Heat transfer and skin friction (or 
viscous) drag are strongly dependent of whether the boundary layer is laminar or turbulent. The 
boundary layer thickness is defined as the distance where the local longitudinal velocity V(y), the 
velocity parallel to the surface of the plate, becomes equal to 99 percent of the freestream velocity V. 

For a flat plate in incompressible flow, the laminar boundary layer thickness Slam can be expressed 
as a function of local Reynolds number and distance from the leading edge of the plate x 

8 L am = 5.0 x/Re x 1/2 


where 


Re x = p V x / (i 


and p and V are freestream density and freestream velocity and p is the absolute viscosity of fluid. 

The skin- friction drag coefficient due to shear stress Cf is equal to iw / q , where xw is the shear 
stress at the body surface and q is the freestream dynamic pressure. In laminar, incompressible flow 
over a flat plate of length L the skin-friction drag coefficient is obtained from the Blasius drag 
equation (see Schlichting, 2000): 


C f lam- L 328 /Re L 1/2 


16 A laminar flow is one in which a fluid flows in parallel layers. In turbulent flow, the fluid velocity vector has unsteady, 
random components added to its mean value. Examples of transition from laminar to turbulent flow are cigarette smoke 
rising in still air and a stream of ink flowing in clear water. In both cases the flow starts out as a uniform, laminar flow. 
As the distance from the point of origin of the flow increases (or as the flow velocity increases) the flow will abruptly 
begin to pass through a transition regime and the smoke or ink will begin to mix, diffuse, and become turbulent. 
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where 


ReL = p V L / p. 


For a flat plate in incompressible flow, the turbulent boundary layer thickness 8turb can also be 
expressed as a function of Re x and x 


Sturb ~ 0.37 x / Re x 1/5 

The skin- friction drag coefficient for a flat plate of length L in turbulent flow is 

Cf turb = 0.074 / ReL 1/5 


If the transition from laminar to turbulent flow is assumed to occur at Rex « 500,000, then the skin 
friction drag coefficient can be estimated from the Prandtl-Schlichting formula for a smooth flat plate 
at zero angle of attack: 


Cf TRANS = 0.455 / (log Re L ) 2 - 58 - 1700 / Re L 


for 5 x 10 5 < Re L < 10 8 

These relations are shown in figure 8 (see Schlichting, 2000 and White, 2005). The skin friction drag 
force can be determined from the equation: 

skin friction drag force = Cf q Sw 
where Sw is the wetted surface area. 

It is very difficult to predict the location and extent of the transition region. Usually a transition line 
location is assumed based on an empirical correlation from experimental data for similar body shapes 
and conditions. The correlation is usually between Reynolds number based on distance from the 
leading edge of a wing or the stagnation point on a blunt body, and the Mach number at the edge of 
the boundary layer. Upstream of the transition line the flow is assumed laminar and downstream the 
flow is assumed turbulent. The location of the transition line is very important for accurate prediction 
of skin friction drag and aerodynamic heating. Anderson (1989) notes that uncertainties in the 
transition region location can have a very significant effect on vehicle weight estimates. In addition 
to surface roughness and freestream turbulence, other factors that can influence boundary layer 
transition are Mach number, Reynolds number, body shape (blunt versus slender) and nose radius, 
wall and stagnation temperatures, ablation (mass removal at the surface), angle of attack, axial 
pressure gradient (decreasing pressure promotes transition), and chemical reactions in the flow 
(Tauber, 1990). 

The variation of skin friction coefficient and Stanton number with Reynolds number is similar to that 
shown in figure 8 for skin friction drag coefficient. There are a variety of methods for predicting skin 
friction in high-speed flow (see Liepmann and Roshko, 1957; Anderson, 1989; and White, 2005). For 
a given Reynolds number, the skin friction coefficient generally decreases as Mach number 
increases. 
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1.6 Flight and Ground Testing 

Flight and ground testing are used to determine vehicle aerodynamic and aerothermodynamic 
characteristics. Flight tests to simulate SRV atmospheric entry conditions are relatively expensive 
and difficult to perform (see sec. 7). One possibility is to launch a model vehicle into the atmosphere 
with a propulsive device. Ground testing is more often employed to simulate SRV flight conditions. 
Many different types of ground test facilities have been used to obtain vehicle pressure distributions 
and aerodynamic force and moment 17 characteristics. Heat transfer measurements and flow 
visualization are possible in conventional ground test facilities. In a ballistic range, data are obtained 
from a model launched from a gun into a quiescent gas (Chapman and Yates, 1998). A shock tube 
can be used to achieve higher stagnation temperatures than a conventional wind tunnel. However, 
shock tube test durations are extremely short and data acquisition is much more difficult. Arc-heated 
(or arc jet) test facilities such as those at the NASA Ames Research Center and U.S. Air Force’s 
Arnold Engineering Development Center (AEDC) in Tullahoma, Tennessee are usually used to 
simulate high- temperature effects. 

The principal flight variables to match for atmospheric entry simulation are: 

Mach number (compressibility effects) 

Reynolds number (viscous effects) 

Total enthalpy (thennal energy) 

Ratio of specific heats (isentropic exponent of appropriate test gas) 

Stagnation pressure (or flight altitude) 

Ratio of wall temperature to stagnation (or freestream) temperature 
Viscous interaction parameter 18 (for some hypervelocity flows) 

Gas chemistry (including dissociation, ionization, and ablation effects) 

Knudsen number (for free molecular or rarefied atmosphere effects) 

Prandtl number 19 (heat conduction) 

It is impossible to achieve complete simulation of entry flight conditions in a single test facility. 
Usually an attempt is made to match the most important variables in one or more facilities. In 
conventional wind-tunnel and ballistic-range testing, flight Mach number and Reynolds number can 
often be matched. Low-density wind tunnels are used to simulate free molecular flows by matching 
the flight Knudsen number. For hypervelocity flows, arc -jet facilities are used to evaluate TPS 
materials. In arc -jet testing an attempt is made to match total enthalpy and stagnation pressure. Note 
that arc -jet tests can not provide useful aerodynamic characteristics. 


17 In aerodynamics as in rigid-body mechanics the moment of an aerodynamic force is the tendency of the force to rotate 
a body about a point or an axis. The moment is the product of the force and a moment arm defined as the perpendicular 
distance from the point or the axis to the line of action of the force. Moments are vector quantities having both 
magnitude and direction. 

18 Viscous interaction parameter is defined in section 3. 

19 The Prandtl number is the ratio of momentum transfer to heat transfer. Prandtl number is defined in section 3. 
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1.7 Computational Fluid Dynamics and the Navier-Stokes Equations 

The Navier-Stokes equations describe the motion of a viscous fluid. Changes in momentum, or 
acceleration, of fluid elements are related to changes in pressure and viscous forces to determine the 
motion and properties of the fluid. The equations can be derived from the basic principles of 
conservation of mass, momentum, and energy. The Navier-Stokes equations are extremely useful in 
the study of fluid dynamic flows; they are higher-order partial differential equations and are very 
difficult to solve. There are limitations and assumptions involved in the derivation; e.g., the fluid 
must be continuous and homogeneous (it does not contain voids or bubbles). In order to obtain useful 
solutions to practical problems it is necessary to use computational fluid dynamics (CFD) techniques. 
For an inviscid (nonviscous) fluid the Navier-Stokes equations reduce to the Euler equations. 20 

Computational fluid dynamics methods are used to obtain detailed flow field characteristics for 
SRVs. The calculated surface pressures and temperatures are used to determine aerodynamic forces 
and heat loads. However, it is absolutely necessary to validate and calibrate CFD methods and 
solutions. The necessary validation and calibration can be accomplished by means of flight or ground 
test. 

Computational methods have been used in fluid dynamics for several decades. Early efforts were for 
inviscid, low-speed and supersonic flows over simple bodies. In the 1980s high-speed computers 
(and a greater understanding of the stability of numerical methods) made possible the solution of 
viscous flows around complicated body shapes (CFD-Online, 2006). 

In CFD the space around the body is divided into small cells to form a mesh or grid and super 
computers are used to iteratively solve the compressible or incompressible, laminar or turbulent fluid 
flow equations. The equations for chemical reactions, species concentration, and heat transfer can 
also be included (Anderson, 1995). CFD solutions are very sensitive to the mesh configuration and 
turbulence model selected. 

Note that there are also numerous computer codes available for heat transfer modeling, structural 
design, trajectory analysis, chemistry and thermodynamic properties. Many of these codes were 
originally developed by NASA and its U.S. space program contractors (Open Channel Software, 
2006). 21 


20 NASA Glenn Research Center has a Learning Technology website that describes the Navier-Stokes equations and 
includes many other subjects in aeronautics (see NASA Glenn Research Center). 

21 The Open Channel Foundation has an agreement with the National Technology Transfer Center to maintain NASA’s 
COSMIC software collection. This collection covers a wide range of disciplines including engineering, chemistry, and 
aerodynamics. 
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2. EARTH ATMOSPHERE AND ACCELERATION DUE TO GRAVITY 


2.1 Properties of Earth’s Atmosphere 

Knowledge of the properties of Earth’s atmosphere is necessary to determine aerodynamic and 
thennal loads during entry. Pressure, temperature, mass density, species concentration, molecular 
weight, viscosity, thennal conductivity and, at higher altitudes, mean free path and collision 
frequency must be known as functions of altitude. Below an altitude of 100 km Earth’s atmosphere is 
primarily nitrogen and oxygen molecules (0 2 dissociation becomes significant at altitudes above 
100 km); see table 3 for 1976 U.S. Standard Atmosphere values (National Oceanic and Atmospheric 
Administration, 1976). 


TABLE 3. PRIMARY CONSTITUENTS OF EARTH’S ATMOSPHERE IN MOLE (VOLUME) 
FRACTIONS, 1976 U.S. STANDARD ATMOSPHERE* 


Altitude, y, km 

n 2 

o 2 

O 

Ar 

0 to 85 

0.7808 

0.2095 

0 

0.0093 

100 

0.7750 

0.1798 

0.0354 

0.0079 

125 

0.7085 

0.0777 

0.2111 

0.0022 

150 

0.6024 

0.0530 

0.3430 

0.0010 


* National Oceanic and Atmospheric Administration, 1976. 


Properties of the atmosphere vary with time of day, day of the year and solar radiation intensity. 
However, these distinctions are usually ignored for hypervelocity Earth entry (Regan, 1984) and a 
standard atmosphere is adopted for trajectory simulation and load calculation. Above an altitude of 
100 km, atmospheric density (and drag force) is so low that its influence on entry vehicle dynamics is 
usually considered negligible. 

Properties of the Earth’s atmosphere at sea level are given in table 4 (National Oceanic and 
Atmospheric Administration, 1976). Properties for the 1976 U.S. Standard Atmosphere can be 
obtained for altitudes up to 150 km from a convenient and easy to use Internet database 
(Squire, 2006). 
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TABLE 4. SEA LEVEL PROPERTIES OF EARTH’S ATMOSPHERE, 
1976 U.S. STANDARD ATMOSPHERE* 


Quantity 

Metric Units 

English Units 

Atmospheric density, p 0 

1.2250 kg/m 3 

0.0765 lb/ft 3 

Atmospheric pressure, p 0 

101,325 Pa 

2116 lb/ft 2 

Atmospheric temperature, T 0 

288.15 K 

518.67 °R 

Mean free path, L 0 

6.6328 x 10" 8 m 

21.76 x 1 0“ 8 ft 

Molecular weight, M 0 

28.9644 kg/kmol 

— 

Speed of sound, Co 

340.294 m/s 

11 16 ft/s 

Thermal conductivity, k 0 

0.025326 W/(m-K) 

0.0146 Btu-ft/ (hr-ft 2 -°R) 

Absolute viscosity, p 0 

1.789 x 10 3 kg/(m-s) 

1.202 x 10 5 lb/(ft-s) 

Acceleration due to gravity, g 0 

9.80665 m/s 2 

32.174 ft/s 2 


* National Oceanic and Atmospheric Administration, 1976. 


For Earth’s atmosphere the mass density of air p and acceleration due to gravity g can be expressed 
as functions of altitude (see secs. 2.1.2 and 2.4). 

Static pressure p can be determined by summing the pressure and gravitational forces acting on a 
fluid element: 


dp/dy = - p g 
P = Po - /p g dy 

where geometric altitude y is measured from the mean Earth radius Ro (see sec. 2.4). 
Note that geopotential altitude h is defined as follows: 

go dh = g dy 


or 

h = 1 [g(y) / go] dy. 


2.1.1 Temperature 

The 1976 Earth model atmosphere consists of constant temperature layers and layers with 
temperature varying linearly with altitude. Between sea level and about 100 km there are four distinct 
layers (see fig. 9). 
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Troposphere (below 1 1 km) the temperature gradient (or lapse rate) is -6.5 K/km. 

Stratosphere (from 1 1 km to about 47 km) the temperature gradient varies from 0 to 2.8 K/km. 

Mesosphere (from about 47 km to 80 km) the temperature gradient varies from 0 to -3.9 K/km. 

Thermosphere (above 80 km) the temperature gradient varies from 0 to 9 K/km. 

The bounding altitudes separating the four regions are referred to as: the Tropopause (11 km); the 
Stratopause (about 47 km); and the Mesopause (80 km). 

2.1.2 Density 

This simple expression is often used to detennine the atmospheric mass density p as a function of 
altitude y 

p(y) = Ke' Py 

where K and P are constants selected to match the preferred model atmosphere; p is called the 
inverse density scale height. 

Allen and Eggers (1958) and Eggers et al. (1958), used values of 0.0034 slugs/ft 3 and 22,000 ft for K 
and 1/p, respectively. Tauber (1998) recommends K = p 0 = 1.225 kg/nr 3 and 1/p = 7200 m. Regan 
and Anandakrishnan (1993) recommend Allen and Eggers' values equal to 1.752 kg/ nr and 6700 m 
for K and 1/p, respectively. Figure 10 shows a comparison between the analytic expression given 
above for p(y) and the 1976 U.S. standard atmosphere. 

2.1.3 Mean Free Path 

The mean free path X is the mean value of the distance traveled by a neutral particle between 
successive collisions with other particles (National Oceanic and Atmospheric Administration, 1976). 
Mean free path is shown as a function of altitude in figure 1 1 . The limits of the continuum and free 
molecular flow regimes are determined by the ratio of the mean free path to a characteristic body 
length. This ratio is called the Knudsen number. 


2.2 Viscosity 

The transport property known as absolute, or dynamic, viscosity p is important in the determination 
of Reynolds numbers and frictional forces. For air at altitudes less than 85 km, absolute viscosity p 
can be determined from the air temperature T by the Sutherland approximation 

p = 1.458 x IQ' 6 T 372 / (T + 1 10.4) kg/m-s 


where T is in Kelvin. 
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2.3 Thermal Conductivity 

The coefficient of thennal conductivity k is an important transport property in heat transfer. For air at 
altitudes less than 85 km, thermal conductivity k can be determined from the air temperature T by a 
formula similar to the Sutherland approximation for viscosity 

k = 2.646 x 1CT 3 T 3/2 / (T + 245.4 x 1CT 12/t ) W/m-K 

where T is in Kelvin (Regan and Anandakrishnan, 1993). 


2.4 Acceleration Due to Gravity 

The acceleration due to gravity g can be expressed as a function of altitude y 

g(y) / go = [Ro / (Ro + y)] 2 = ( l + y / Rof 2 
~ [ 1 _ 2 (y / R 0 ) + 3 (y / R 0 ) 2 + . . . ] 

where 

go = g(0) = 9.80665 m/s 2 


and 


Ro = mean or effective radius of the Earth = 6,356,766 m. 


The above approximation for g(y) / go is derived from the binomial theorem for y / Ro « 1 (see 
Zwillinger, 2002). Note that the Earth's mean equatorial radius is 6,378,140 m (approximately 21 km 
larger than its polar radius). The Earth's equatorial bulge and its nonuniform mass density result in 
deviations from the inverse- square approximation for spherical bodies. These deviations are usually 
accounted for in trajectory simulation codes. The Earth’s oblateness is accounted for by including the 
2nd zonal harmonic term h in the gravitational model. 


The International Gravity Formula gives g at sea level as a function of latitude X 
go = 9.780495 [1 + 0.0052892 sin 2 (k) - 0.0000073 sin 2 (2 X)] 


For X = 45 deg this fonnula gives: go = 9.806289 m/s 2 (Glover, 2006). 
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3. AERODYNAMICS AND HEAT TRANSFER 


3.1 Basic Aerodynamic Concepts 

The importance of dimensionless parameters, such as Mach number, Reynolds number and Knudsen 
number, in the study of atmospheric entry was mentioned in section 1 . Dimensionless force and 
moment coefficients are equally important in trajectory computation. If a vehicle’s surface pressure 
distribution is known, the aerodynamic forces and moments can be determined by integrating the 
pressure over the surface area. Often it is possible to test a scale model of a vehicle in a wind tunnel 
and measure the forces and moments on the model. If the test Mach number and the test Reynolds 
number match those in flight, then the forces and moments can be scaled to flight conditions. 22 (If 
free molecular flow will be encountered in flight, then ground-test Knudsen number should match the 
flight Knudsen number.) 

An aerodynamic force coefficient is defined as 

Force coefficient = Force / (dynamic pressure x reference area) 

A moment coefficient is defined as 

Moment coefficient = Moment / (dynamic pressure x reference area x reference length) 


where 


dynamic pressure = q = 1/2 p V 2 
or for a thermally perfect gas (not valid for M » 1) 

q = 1/2 y p M 2 

p = freestream mass density of fluid 
V = freestream fluid velocity 
M= freestream Mach number 
P = freestream static pressure 
y = isentropic exponent. 


22 Force Coefficient = Force FUG HT / ^flight x SfLIGHt) - Force WIND TUNNEL / (eiwiND TUNNEL x S WIND TUNNEl) 
or ForcepLiGHT = ForcewiND tunnel (^flight / c[wind tunnel) (Sflight / Swind tunnel) = Force Coefficient x q FL iGHT x Sflight- 
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In a wind tunnel test a calibrated strain gauge balance is used to measure forces and moments. The 
axial force, normal force, and pitching moment are of principal concern: 

Axial force = A = C A q S 

Normal force = N = Cn q S 

Pitching moment = M = C m q S L 


or 


Axial force coefficient = C A = A / q S 
Normal force coefficient = Cn = N / q S 
Pitching moment coefficient = C m = M/ qSL 

where 

S = reference area (cross-sectional or planform area) 

L = reference length. 

Once C A and Cn are known it is a simple matter to determine the lift and drag coefficients and the lift 
and drag forces 


Cl = Cn cos a - C a sin a 
Cd = Cn sin a + C A cos a 

where a = angle of attack (the angle between the vehicle axis and the velocity vector) 

Lift force = Cl q S 
Drag force = Cd q S. 

Figure 12 shows how axial and normal forces are resolved into lift and drag by vector addition. Note 
that the vector sum R = L + D = N + A. 

For some body shapes, analytical and/or experimental lift, drag and moment coefficient data are 
readily available for the entire Mach number range from subsonic to hypersonic. These data can be 
found in handbooks and textbooks; see Hoerner (1965 and 1985), McConnick (1995), and Avallone 
and Baumeister (1996). 
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3.2 Aerodynamic Drag 

The aerodynamic force on a vehicle acting directly opposite to the vehicle’s velocity vector is called 
drag. In the absence of thrust, atmospheric drag results in vehicle deceleration. Various components 
of drag force are described below. Drag estimation can be a very involved process. 

Pressure or form drag (usually called wave drag in supersonic or hypersonic flow) is the result of 
pressure forces acting on the vehicle surface. In hypersonic flow, the Newtonian approximation 
(described below) usually provides an adequate estimate of the pressure drag. 

Skin friction drag is the result of viscous shear stresses acting on the vehicle surface. The skin 
friction drag coefficient Cf is based on the vehicle surface (or “wetted”) area Sw and not on a cross- 
sectional area. 23 The skin friction drag coefficient increases when the flow in the boundary layer 
changes from laminar to turbulent. Modem fighter aircraft have skin friction drag coefficients around 
0.003 (McCormick, 1995). The skin friction drag coefficient for a flat plate is shown as a function of 
Reynolds number in Figure 8. Note that at Rex = 1,000,000, the incompressible turbulent skin 
friction drag coefficient is 3 times greater than the incompressible laminar skin friction drag 
coefficient. 24 

Base drag is the drag produced by the pressure acting on the base of the vehicle and can be negative 
(acting in the upstream direction). Base drag is a function of the body and afterbody shapes, Mach 
number, and Reynolds number. 25 

Another way to look at drag is to consider it to be composed of a component independent of angle of 
attack (zero-lift drag) and a component that is a function of angle of attack (called induced drag). For 
wing sections, induced drag usually increases as the square of the lift. 

In subsonic flow, drag is relatively unaffected by Mach number. As Mach number approaches one, 
shock waves begin to fonn and the drag coefficient begins to increase significantly (transonic flow 
regime). In supersonic flow, the drag coefficient may continue to increase, remain constant, or 
decrease depending on the body shape (McConnick, 1995). Typical values of drag coefficient are 
shown in table 5. 

The drag on two-dimensional bodies (e.g., infinite cylinders with axis normal to the freestream 
direction and wedges of infinite width) is usually greater than the drag of corresponding three- 
dimensional bodies (spheres and cones). 


23 When adding skin friction and pressure drag coefficients the skin friction coefficient must first be multiplied by the 
ratio of the surface area to the cross-sectional area Sw/S. 

24 Convective heat transfer varies directly as skin friction and follows a similar trend as Reynolds number is increased. 

25 The drag force measured in a wind tunnel may not include base drag if the wind-tunnel model support (usually called a 
“sting”) is attached at the base of the model. 
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TABLE 5. TYPICAL VALUES OF DRAG COEFFICIENT* 


Body Shape 

Incompressible 
p = constant 
M«0 

Hypersonic 
M » 1 

Flat Plate 

Surface parallel to velocity vector, 
a = 0 (skin friction drag; see fig. 8) 

1.33/(Re L ) 1/2 

(laminar flow) 

0.074/(Re L ) 1/5 

(turbulent flow) 

— 

Surface perpendicular to velocity 
vector, a = 90 deg 

1.2 

(for F « width) 

p 

<N 

1 

00 

Circular Disc 

1.2 

p 

<N 

1 

OC 

Circular Cylinder 

0.6 (for F « diameter) 

1.3 

(axis perpendicular to velocity vector) 

1000 <Re D < 100,000 

1 .2 (for L » diameter) 

— 

Sphere 

1000 <Re D < 300,000 

0.4 

1.0 

Re D * 500,000 

0.2 

— 

Cone 

Semivertex angle 0c = 15 deg 

0.34 

0. 1 34 a 

Semivertex angle 0 C = 30 deg 

0.51 

0.500 a 


*Avallone and Baumeister, 1996; McCormick, 1995; Hughes and Brighton, 1999; and White, 2005. 
a Inviscid hypersonic flow (does not include skin friction drag), perfect gas, y =1.4 (see table 6). 


3.3 Newtonian Impact Theory for Inviscid Hypersonic Flow 

A law proposed by Isaac Newton more than three centuries ago can be used to predict lift and drag in 

inviscid (nonviscous) hypersonic flow (Anderson, 1989). Newton’s sine-squared law states that the 

2 2 

force acting on an inclined body surface is equal to p V“ S sin" 0 

F = p V 2 S sin 2 0 = A p S 


where 

0 = angle between tangent to body surface and freestream velocity vector 
A p = surface pressure - freestream static pressure 
S = surface area 
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or expressed in terms of a pressure coefficient C p 


C p = Ap/q = 2 sin 2 0 


where 

q = 1/2 p V 2 . 

This result is of no value in subsonic or supersonic flows and is not applicable to downstream facing 
surfaces (where 0 < 0 and C p is assumed equal to zero). Inviscid hypersonic aerodynamic force 
coefficients can be obtained by integrating the Newtonian pressure coefficient over the body surface. 

Newtonian lift and drag coefficients Cl and Cd for a flat plate at angle of attack a are 

Cl = 2 sin 2 a cos a 

C D = 2 sin 3 a 

L/D = cot a 

These quantities are shown as functions of angle of attack in figure 13. Note that the lift coefficient 
has a maximum value of 0.770 at a = 54.7 deg. The value of L/D increases as a decreases. However, 
L/D will go to zero as a and L go to zero because D has a finite value at a = 0 due to skin friction 
drag. 

In inviscid hypersonic flow the Newtonian approximation gives the following values for the wave or 
pressure drag coefficient: 

For a circular cylinder with its axis perpendicular to the freestream direction 

C D = 4/3 


For a sphere 


C D =1 


For a cone with semivertex angle 0c 

Cd = 2 shr 0c 

A modification to Newtonian theory proposed by Lester Lees gives more accurate results for blunt 
bodies. Lees' approximation replaces the coefficient 2 in the sine-squared law with the pressure 
coefficient at the stagnation point behind the nonnal shock wave C p max (see Anderson, 1989) 

Cp C p max sin 0 

For a perfect gas, C p max can be expressed as a function of freestream Mach number M and the 
isentropic exponent y (Ames Research Staff, 1953). 

As M goes to oo C p max goes to 1 .839 (for y = 1.4) 
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The following equation gives the Newtonian drag coefficient Cd for a sphere-cone at zero angle of 
attack with nose radius r N , base radius re, and cone-section semivertex angle 9c : 


C D = 2 sin 2 9c + (r N / r B ) 2 cos 4 9c 

Cd is shown as a function of bluntness ratio r\j / r B and cone-section semivertex angle 9c in figure 14. 
(When lift is zero, the drag coefficient is referred to as the zero-lift drag coefficient; e.g., axially- 
symmetric bodies have zero lift at zero angle of attack). 

The lift and drag coefficients presented above are for forces due to surface pressure in inviscid 
hypersonic flow. In supersonic and hypersonic flow the pressure drag is called wave drag because the 
surface pressure increase is due to flow compression as the flow passes through the shock wave. 

Hypersonic lift, drag and pressure coefficients for many body shapes can be readily obtained from 
simple computer algorithms (Regan and Anandakrishnan, 1993). Note that in this section viscous 
effects are neglected and only inviscid pressure forces are considered. 


3.4 Other Approximate Methods for Inviscid Hypersonic Flow 

The static pressure at any point on the surface of a body in hypersonic flow can often be determined 
from the slope of the body surface at that point. These approximations are referred to as local surface 
inclination methods (see Anderson, 1989). 

3.4.1 Tangent-Wedge Method 

The tangent-wedge method can be applied to sharp-nosed, two-dimensional bodies with attached 
shock waves. The local static pressure at a point on the body surface is assumed equal to the static 
pressure behind an oblique shock wave with deflection angle equal to the slope of the body surface at 
the point in question; i.e., the equivalent wedge has deflection angle equal to the local body slope. 

The local static pressure is a function of freestream Mach number M, the local deflection angle 9 and 
the isentropic exponent y. Lift and drag can be obtained by integrating the pressure over the body 
surface. For SRVs two-dimensional body sections are likely to occur only as appendages to the basic 
body (e.g., wings, fins, external panels, etc.). Flow properties downstream of oblique shock waves 
are given in NACA Report 1 135 (Ames Research Staff, 1953); an Internet calculator is also available 
(Devenport, 2996). 

3.4.2 Tangent-Cone Method 

The tangent-cone method for sharp-nosed, three-dimensional (or axially-symmetric) bodies is similar 
to the tangent-wedge method for two-dimensional bodies. The tangent-cone method also applies only 
to sharp-nosed bodies with attached shock waves. Static pressure on the body surface is obtained 
from supersonic conical flow solutions (see Sims, 1964). This reference tabulates the flow properties 
between the conical shock wave and the cone surface. The tabulated values were obtained from 
numerical solutions of ordinary differential equations for conical flow with the ratio of specific heats, 
freestream Mach number, and cone semivertex angle as parameters. Both tangent-cone and tangent- 
wedge methods give reasonably good results for lift, drag and pressure on sharp-nosed bodies. 
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3.4.3 Shock-Expansion Method 


For a body with an attached, or a detached shock wave the flow is first compressed as it passes 
through the shock wave and it then expands as it flows along the body surface (assuming convex 
body curvature). A uniform, two-dimensional, supersonic flow over a convex curve is known as a 
Prandtl-Meyer expansion fan. For a Prandtl-Meyer expansion, the local Mach number is a function of 
the initial Mach number M, the turning or deflection angle and the isentropic exponent y. The local 
static pressure and pressure coefficient can be determined from the local Mach number and the 
isentropic flow relations; see NACA Report 1135. Lift and drag are obtained by integrating the 
pressure over the body surface. Although the Prandtl-Meyer expansion is intrinsically two- 
dimensional it provides reasonably good results for axially- symmetric bodies in hypersonic flow. An 
Internet calculator is available for obtaining Prandtl-Meyer expansion flow properties (Delft 
Technical University, 2006). 


3.5 Shock Wave Detachment 

Sharp-nosed bodies usually have attached shock waves in supersonic and hypersonic flow (oblique 
shock waves for two-dimensional wedge-shaped bodies and conical shock waves for conical bodies). 
If the freestream Mach number is decreased or if the apex angle of a pointed body is increased the 
shock wave will detach and the flow will be similar to that for the blunt body shown in figure 1 . 
Figure 15 shows sharp-nosed bodies with attached and detached shock waves. (Blunt bodies have 
detached shock waves in supersonic and hypersonic flow.) Minimum Mach numbers and maximum 
semivertex angles for attached shock waves for cones and wedges in inviscid, perfect gas flow are 
given in NACA Report 1 135. The conical flow charts in NACA Report 1 135 are for y = 1.405 
(Kopal, 1947). Conical flow solutions for y = 1.4 are given in NASA SP-3004 (Sims, 1964). NASA 
SP-3007 contains solutions for cones at small angles of attack (Sims, 1964). 


3.6 Hypersonic Mach Number Independence 

Analytical solutions and experimental results show that for very high freestream Mach numbers, 
many flow characteristics (including lift, drag, pressure and moment coefficients) are relatively 
insensitive to changes in Mach number. This phenomenon is observed in experimental data at 
hypersonic Mach numbers for both sharp and blunt, two- and three-dimensional bodies. For sharp 
circular cones in inviscid, perfect gas flow, table 6 shows that the zero-lift drag coefficient 
approaches the Newtonian values as the Mach number increases and the change in drag coefficient is 
very small for freestream Mach numbers greater than ten. 
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TABLE 6. ZERO-LIFT DRAG COEFFICIENT C Dq FOR SHARP, CIRCULAR CONES, 

INVISCID, PERFECT GAS WITH y = 1.4* 


Freestream 
Mach number, 
M 

Cone semivertex angle, 0c, degrees 

10 

20 

30 

1.5 

0.124 

0.387 

b 

2.5 

0.094 

0.299 

0.598 

3.5 

0.083 

0.275 

0.561 

5 

0.075 

0.260 

0.542 

7 

0.070 

0.253 

0.533 

10 

0.067 

0.249 

0.529 

15 

0.065 

0.247 

0.526 

20 

0.064 

0.246 

0.525 

Newtonian 3 

0.060 

0.234 

0.500 


* NASA SP-3004 (Sims, 1964). 

a The Newtonian zero-lift drag coefficient is equal to 2 sin 2 9 C . (For right circular cones at zero angle of attack the 
zero-lift drag coefficient is equal to the pressure coefficient.) 
b Shock wave is detached. 


3.7 Hypervelocity and High Temperature Effects 
3.7.1 Wall Catalysis 

In chemically reacting flows over solid surfaces some surface materials promote recombination of 
atoms. Surfaces that promote complete recombination of all atoms are called catalytic. (If no 
recombination occurs at the surface, the surface is referred to as noncatalytic.) Wall catalysis is a 
concern because recombination reactions can increase the local heat flux by a factor of 2 or more. It 
is also possible to have partially catalytic surfaces that are very difficult to model because of the 
complexity of nonequilibrium chemistry. Generally a vehicle designer will attempt to select a 
noncatalytic surface (a surface which does not promote recombination reactions). 

Lightweight ceramic tiles with reaction cured glass (RCG) coatings are characterized as noncatalytic 
TPS materials. Other coatings, such as toughened uni-piece fibrous insulation (TUFI) have been 
developed to provide low catalycity and increased strength. Metallic surfaces, on the other hand, are 
highly catalytic and promote recombination reactions. Metallic TPS materials have been proposed 
(e.g., for reusable launch vehicles) but suffer from the problem of high catalycity as well as high 
thennal expansion (this makes noncatalytic surface attachment or coating difficult). Partially catalytic 
TPS materials include reinforced carbon carbon (RCC) and silicon carbide materials. 
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When using CFD codes, surface catalysis is an important boundary condition and the mass fraction 
of each chemical species is a dependent variable. For a fully catalytic wall, recombination reactions 
occur at an infinite rate. 

3.7.2 Viscous Interaction 

In classical boundary layer theory, it is assumed that viscous effects are confined to a small region 
adjacent to the surface of the body. In hypersonic flow, boundary layer thickness increases as the 
square of the Mach number. Thus for high Mach number flows, the classical small boundary layer 
thickness assumption (8 « L) is no longer valid and there is significant interaction between the 
viscous (boundary layer) and the nonviscous regions of the flow field. This is referred to as the 
viscous merged layer flow regime and it occurs at higher altitudes and velocities than continuum 
flow. 

Hypersonic vehicle surface pressures can be correlated using a viscous interaction parameter % 

X = M 3 (C/Re) 1/2 


where 

C = Chapman-Rubesin constant = pw p w /p e p e 
Subscript W indicates wall or surface value 
Subscript e indicates value at outer edge of boundary layer. 


Near the body leading edge or stagnation point (x > 3) there is a strong viscous interaction. Viscous 
interaction is weaker for smaller values of %. Pressure and force coefficients (as well as skin friction 
and heat transfer coefficients) can be correlated with a modified viscous interaction parameter V 

V = M (C/Re) 1/2 


For hypersonic vehicles the maximum L/D decreases as the modified viscous interaction parameter 
increases. This is because the change in lift (with increase in V) is small compared with the increase 
in pressure and skin friction drag. 

Another type of viscous interaction occurs when a shock wave intersects or impinges on a boundary 
layer; this is called shock-boundary layer interaction. The incident shock can cause local flow 
separation. The separated boundary layer may induce a second shock wave. When the boundary layer 
reattaches to the surface a third, or reattachment, shock wave may be formed (Anderson, 1989). In 
hypersonic flows, shock-boundary layer interaction can cause severe local aerodynamic heating 
problems. 

Local heat fluxes can increase dramatically due to shock wave impingement. Hypersonic pressure 
and heat transfer effects in complex flow geometries can be extremely difficult to predict and 
experimental verification is normally required. NASA’s Mach 6.7 experimental X-15 aircraft 
experienced severe heating caused by shock wave impingement and resulted in pylon failure and 
damage to the aircraft. 
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3.8 Heat Transfer 


Structural materials, and materials used for thermal protection, have a maximum allowable 
temperature and can absorb a limited quantity of heat. The vehicle must also be designed to 
withstand the thennal stresses induced by temperature gradients as heat flux can vary along the 
vehicle surface. Structural strength decreases as temperature and heat flux increase. Excessively high 
heat fluxes can cause spalling (loss of surface material) or structural failure due to thennal stress. 

Heat is transferred from the high-temperature shock layer to the vehicle by forced convection and by 
radiation. Some heat may be re -radiated to space. Heat is transferred by conduction through the 
structure to the vehicle interior. A circulating liquid or transpiration cooling system (injecting cold 
gas into the boundary layer) can dissipate some of this heat. An active cooling system may not be 
practical for an SRV due to added weight and complexity. If an ablative heat shield is employed, 
other modes of heat and mass transfer must be considered. (Ablation is discussed in sec. 6.) 

3.8.1 Convective Heat Transfer 

Many of the approximations for convective heat transfer in high-temperature gasses are based on the 
work of Fay and Riddell, 1958. The convective heat flux at the stagnation point qcoNvt can be 
expressed as a function of freestream velocity V, freestream mass density p, and the body radius, or 
radius of curvature, r N (Tauber, 1989) 

qcoNvt * 1.83 x 10 4 (p/r N ) 1/2 V 3 [1 - (h w / h t )] W/m 2 


where 

h t = stagnation enthalpy = c p T + 1/2 V (for a perfect gas) 
hw = specific enthalpy at the vehicle surface 
p is in kg/m , r N is in m, V is in m/s, and h is in J/kg 


Enthalpy is a function of wall temperature Tw and stagnation pressure p t . Total enthalpy, or heat 
content, is the sum of internal energy and the product of volume multiplied by the pressure. Specific 
enthalpy, or enthalpy per unit mass, h is the sum of internal energy per unit mass and the ratio of the 
pressure divided by density (Liepmann and Roshko, 1957). Internal energy is the kinetic energy due 
to molecular motion (translational, rotational, and vibrational) and potential energy associated with 
the vibrational and electric energy of atoms within molecules. The change in internal energy is equal 
to the sum of the heat input to a system and the work done on or by the system. 

The Tauber approximation given above is much simpler and usually more conservative than the more 
accurate Fay-Riddell solution. 
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For h w « h t , (cold- wall approximation) the maximum stagnation point convective heat flux 
flcoNv t max for “steep” (flightpath angle Ye < -10 deg), nonlifting entry is (Tauber, 1998) 


flcoNv t max ~ 7.5 x 10 7 Ve 3 (BP sin |ye| / tn) 1 ^ W/m 

where the ballistic parameter 26 BP is equal to m/Co S in kg/m 2 and Ye is the flightpath angle at entry 
(see fig. 16). 

Note that cold- wall heat flux is not a realistic prediction of actual flight conditions. Cold-wall heat 
fluxes should not be used for design purposes even though they may be useful for comparing designs. 

For “steep,” nonlifting trajectories, peak stagnation point convective heat flux occurs at altitude 
y = (1 / p) In [-p 0 / (P BP sin Ye)]; po and p are defined in section 2. 

Note that for a given ballistic parameter and a given nose radius, the peak stagnation point heat flux 
decreases as the flightpath angle increases (negative flightpath angles point downward). Peak 
stagnation point heat flux is inversely proportional to the square root of nose radius and directly 
proportional to entry velocity cubed and the square root of ballistic parameter. 

For h w « h t , the stagnation point convective total heat load Qconv t for “steep,” nonlifting entry is 
the integral of the heat flux 


Qconv t * 0.028 V E 2 BP 1/2 /(r N sin | Y e |) 1/2 J/m 2 

Stagnation point convective total heat load is inversely proportional to the square root of nose radius 
and directly proportional to entry velocity squared and the square root of ballistic parameter. For a 
given ballistic parameter and a given nose radius, stagnation point convective total heat load 
decreases as the flightpath angle becomes steeper. (Flight times decrease as the trajectory becomes 
steeper.) 

A “normalized” cold-wall peak stagnation point convective heat flux 27 qcoNv t max (r N /BP) 1/2 is given 
in table 7 and figure 17 as a function of entry velocity for various flightpath entry angles. 

A “normalized” cold- wall stagnation point total heat load Qconv t On/BP) “ is given in table 8 and 
figure 18 as a function of entry velocity for various flightpath entry angles. 


26 In ballistics the quantity W/CdA is referred to as the ballistic coefficient and is usually given in English units, pounds 
per square foot. In this document the ballistic parameter BP is defined as m/C D S and is given in metric units, kg per 
square meter. (The reference areas S and A are equivalent.) 

2 A normalized quantity can be defined as a function or a variable that is multiplied or divided by one or more basic 
parameters. The normalized quantity is used here to scale the physical quantity (heat flux or total heat load) for various 
parameters such as vehicle mass, size, and shape (drag coefficient) and various initial conditions (entry velocity and 
flightpath angle at entry). 
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TABLE 7. “NORMALIZED” COLD-WALL PEAK STAGNATION POINT CONVECTIVE HEAT 
FLUX, qcoNvtMAx (r N /BP) 1/2 IN kW/(kg-m) 1/2 FOR “STEEP” NONLIFTING ENTRY* 


Entry 

Velocity, V E 
km/s 

Flightpath Angle at Entry, Ye, deg 

-10 

-15 

-25 

-40 

-65 

-90 

11 

416 

508 

649 

800 

950 

998 

12 

540 

659 

843 

1039 

1234 

1246 

13 

687 

838 

1071 

1321 

1569 

1648 

14 

858 

1047 

1338 

1650 

1959 

2058 

15 

1055 

1288 

1646 

2029 

2410 

2531 


* h w « h, (cold- wall approximation); 1/p = 7200 m (see Tauber, 1998). 


Example : 

If Ve = 15 km/s and Ye = —25 deg, then qcoNv t max (r N /BP) 1/2 = 1 646 kW/(kg-m) m . 

If r N = 0.1 m and BP = 10 kg/m 2 , then qcoNv t max = 1646 (BP/r N ) 12 = 1646 (10/0. 1) 1/2 = 
16,460 kW/m 2 = 1646 W/crn 2 . 


TABLE 8. “NORMALIZED” COLD-WALL STAGNATION POINT CONVECTIVE TOTAL 
HEAT LOAD Qconv t(r N /BP) 1/2 IN kJ/(kg-m) 1/2 FOR “STEEP” NONLIFTING ENTRY* 


Entry Flightpath Angle at Entry, Ye, deg 

Velocity, Ve 

n m/c -10 -15 -25 -40 -65 -90 


11 

i 

8130 

6660 

5210 

4230 

3560 

3390 

12 

9680 

7920 

6200 

5030 

4240 

4030 

13 

11,360 

9300 

7280 

5900 

4970 

4730 

14 

13,170 

10,790 

8440 

6840 

5760 

5490 

15 

15,120 

12,380 

9690 

7860 

6620 

6300 


* h w « h, (cold- wall approximation); 1/p = 7200 m (see Tauber, 1998). 


Example : 

If Ve = 15 km/s and Ye = -25 deg, then Qconv tOxf/BP) 1/2 = 9690 kJ/(kg-tn) 1/2 . 

If r N = 0.1 m and BP = 10 kg/m 2 , then Q C0 Nvt= 9690 (BP/r N ) 1/2 = 9690 (10/0.1) 172 = 96,900 kJ/m 2 or 
9690 J/ctn 2 . 
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3.8.2 Radiative Heat Transfer 


At Mars SRV Earth entry speeds the entry vehicle will be subjected to significant amounts of 
radiative heating from the high-temperature shock layer gasses; at these entry speeds the radiative 
heating can exceed the convective heating. High-speed radiative heat transfer is a complex, nonlinear 
process and numerical solutions are very difficult to obtain. A correlation can be used to predict the 
stagnation point radiative heat flux as a function of r N , V and p (Tauber and Sutton, 1991): 

qRADt = Cr N A p B f(V) 


where 

A is a function of p and V 
B and C are constants 

f(V) is a tabulated function for Earth’s atmosphere 
and r N , V and p are as previously defined. 

Note that the stagnation point radiative heat flux increases with increased nose radius. The opposite is 
true for stagnation point convective heat flux which decreases as nose radius is increased. 
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4. DYNAMICS AND FLIGHT MECHANICS 


4.1 Shuttle Orbiter Re-entry 

A sample return vehicle will enter the Earth's atmosphere at a much higher velocity than the Space 
Shuttle Orbiter. - However, the Orbiter is a familiar vehicle and it is useful to examine its re-entry 
from low Earth orbit (LEO) because it exhibits many characteristics similar to an SRV. 

The Space Shuttle Orbiter orbits the Earth at an altitude of 320 km (200 mi) and a velocity of 
7.7 km/s. This is in the free molecular flow regime where the molecular mean free path is large 
compared to a characteristic vehicle dimension. (The mean free path is defined as the average 
distance that gas molecules travel between successive collisions with one another; see section 2.1.3.) 
About 60 minutes before touchdown a de-orbit burn occurs. The re-entry event sequence is shown in 
table 9. The vehicle is turned tail forward and the Orbital Maneuvering System (OMS) engines are 
fired for 2 to 3 minutes. At about 160-km altitude, the vehicle leaves the free molecular flow regime. 
The vehicle continues to descend at hypervelocity speed and reaches the “entry interface altitude” 
(122 km) about 30 minutes after initiation of the de-orbit bum. About 5 minutes later the 
communications blackout begins. This is caused by ionized particles enveloping the vehicle. 
Transition from a laminar to a turbulent boundary layer occurs as Reynolds number increases above 
1,000,000. (This value of transition Reynolds number is appropriate for smooth bodies with convex 
surfaces.) Maximum heating (nose and leading edge temperatures about 1500 °C) occurs about 
20 minutes prior to touchdown. The flow around the vehicle is not in chemical or thermodynamic 
equilibrium and oxygen and nitrogen atoms, NO ions, and electrons are present. Continuum flow 
begins at an altitude of approximately 85 km and a Mach number around 25. Roll reversal (S-turn) 
maneuvers begin 16 minutes prior to touchdown. The communications blackout ends about 
12 minutes prior to touchdown. Other high-temperature effects (e.g., nonequilibrium, molecular 
dissociation and vibrational excitation) become insignificant as the vehicle continues to decelerate. 
Braking begins 10 minutes prior to touchdown. At 6 minutes prior to touchdown the flow is no 
longer hypersonic (M « 3). At 3 minutes prior to touchdown the vehicle is traveling at sonic velocity 
(295 m/s at an altitude of 15 km); this is the brief transonic flow regime. At 2 minutes prior to 
touchdown the flow becomes subsonic (M < 1) and the vehicle is on a 22 deg glide slope. At 
30 seconds prior to touchdown a flare maneuver begins and the glide slope is reduced to 1.5 deg. 
Landing speed is approximately 94 m/s. After the Orbiter has landed, cooling vans are dispatched to 
remove the heat accumulated in the TPS tiles during re-entry. 

The flow regimes for a typical Space Shuttle Orbiter re-entry and landing are shown in figure 19. 
Note that the values and flow regime boundaries presented above and in figure 19 apply only to 
Orbiter-size vehicles and trajectories. Smaller vehicles would encounter free molecular flow at lower 
altitudes. 

The important high velocity and high temperature effects can be summarized as follows (Hansen and 
Heims, 1958; Howe, 1990; and Anderson, 1989): 


28 The Space Shuttle Orbiter is also referred to as the Space Transportation System (STS) Orbiter. 
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The gas is calorically perfect at velocity < 1 km/s or temperature < 800 K 
Vibrational excitation 29 occurs at velocity > 1 km/s (temperature > 800 K) 

Oxygen dissociation occurs at velocity > 2.5 km/s (temperature > 2500 K) 

Nitrogen dissociation and NO ionization occur at velocity > 5 km/s (temperature > 4000 K) 
Ionization of atomic nitrogen and oxygen occurs at velocity >10 km/s (temperature > 9000 K) 

TABLE 9. SPACE SHUTTLE ORBITER RE-ENTRY EVENTS AND FLOW REGIMES 


Re-entry Event and Flow Field Characteristics 

Time before 
Touchdown, 
min 

Velocity, 

km/s 

Altitude, 

km 

On-orbit, free molecular flow (hypervelocity) 

— 

7.7 

320 

De-orbit burn 

60 

7.7 

282 

Pitch up to « 30 deg 

52 

7.6 

— 

Begin transition from free molecular to 
continuum flow 

— 

— 

* 160 

Entry interface, atmospheric effects begin 

30 

7.6 

122 

Blackout begins/NO ionization 

25 

7.4 

95 

Begin continuum flow regime, hypersonic flow 
(M « 25) 

— 

— 

« 85 

Maximum heating (boundary layer 
transition/ maximum dec deration) 3 

20 

6.7 

70 

Dissociation of NO ends, roll reversal begins 

16 

~ 5 

— 

Blackout ends 

12 

3.6 

55 

O 2 dissociation ends, braking begins 

10 

« 2.5 

— 

Vibration excitation ends, supersonic flow (M « 3) 

6 

0.99 

27 

Transonic flow (M « 1) 

3 

0.30 

15 

Subsonic flow (M < 1) 

— 

— 

— 

Hold 22 deg glide slope 

2 

0.19 

4.0 

Flare begins (1.5 deg glide slope) 

0.5 

0.16 

0.6 

Touchdown (or recovery) 

0 

0.09 

0 


a The values shown are for maximum heating. Boundary layer transition generally occurs earlier than maximum 
heating; peak deceleration generally occurs later than maximum heating. 


29 


The vibrational energy of the molecules becomes excited and the specific heats and the isentropic exponent are no 
longer constant. 
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4.2 Trajectory Parameters and Vehicle Design Considerations 

In 1953 Allen and Eggers 30 concluded in their study of ballistic missiles entering the Earth’s 
atmosphere at high speeds (Allen and Eggers, 1958): “Blunt shapes appear superior to slender shapes 
from the standpoint of having lower maximum convective heat-transfer rates in the region of the 
nose.” This became the basis for NASA’s early spacecraft. The Mercury, Gemini and Apollo re-entry 
spacecraft had blunt, ablative heat shields. 

Other conclusions about hypervelocity re-entry reached by Allen and Eggers include: a) high-drag 
(sphere-like) bodies should be selected to minimize convective heat transfer (large nose radius); 
b) vehicle weight, surface area, and skin friction drag should be minimized; c) vehicle surface should 
be stiff, strong and curved; d) aerodynamic force sensitivity to angle of attack should be minimized; 
and e) the body should be rotated in flight to equalize heat transfer to all surface elements, reduce 
temperature gradients, and achieve a more unifonn surface temperature distribution. 

Allen and Eggers recognized that high drag and low terminal speeds could result in excessive wind 
drift and high peak deceleration. They noted that these disadvantages could be alleviated with a 
variable geometry extendible skirt at the base of the vehicle to modulate drag. 

SRVs and SRV trajectories must be designed to ensure survival of the vehicle payload after 
atmospheric entry. Special consideration must be given to maximum heat flux, total heat load (the 
integral of heat flux), and maximum deceleration. Peak structural loads generally occur at peak 
dynamic pressure or when the product of dynamic pressure and angle of attack reaches a maximum. 
These are some of the aerodynamic and aerothennodynamic factors that influence SRV design. 


4.3 Equations of Motion 

The equations of motion for a vehicle entering a planetary atmosphere can be written in inertial 
coordinates for a rotating planet. It is much simpler to assume planar motion over a nonrotating 
planet. For Earth the inertial velocity differs from the relative velocity by a westward-pointing vector 
of magnitude 462.8 cos X m/s at sea level, where X is the Earth latitude. 31 

If vehicle mass m is assumed constant (no thrust or mass loss due to ablation), then (Ashley, 1992) 
m dV/dt = - D - in g sin y 
m V dy/dt = (1/r) m V cos y + L- mg cos y 
dr/dt = dy/dt = V sin y 

where r is the distance from the center of the planet to the vehicle center of mass (see fig. 16). 


30 The original report was declassified in 1957 and published in 1958 as NACA Report 1381; see NACA-TR-1381 in the 
NACA report online database (National Advisory Committee for Aeronautics, 2006). 

31 Ashley (1992) claims errors due to neglecting Earth rotation are usually much less than 10 percent. 
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In terms of the ballistic parameter BP and the lift-drag ratio L/D: 
dV/dt = -p V 2 / (2 BP) - g sin y 
dy/dt = (L/D) p V / (2 BP) - (g - V 2 / r) cos y / V 
dy/dt = V sin y 

where (see table 4 and section 2.4) 

BP = m/C D S 
g = go (1 + y/RoT 2 
go = 9.80665 m/s 2 
r = R 0 + y 
Ro = 6,356,766 m 
p ~ p 0 e ' Py 
po = 1.2250 kg/m 3 
p = 1/7200 m 1 

Approximate analytical solutions of these equations can be readily obtained (Allen and Eggers, 1958 
and Chapman, 1959 among others). These equations can also be solved by simple numerical 
integration for given initial conditions with BP and L/D assumed constant. Reasonably good results 
can be obtained by using a fourth-order Runge-Kutta method (Scarborough, 1966) with sufficiently 
small time steps (1 second or less). Runge-Kutta-Fehlberg methods vary the size of the time step 
depending on the behavior of the solution; step size is reduced if the solution varies rapidly 
(Fehlberg, 1969 and Forsythe et ah, 1977). 32 


4.4 Maximum Deceleration for Steep Nonlifting Entry 


For “steep” (y E c -10 deg), nonlifting entry (L/D = 0), with constant drag (BP = const), constant 
gravitational acceleration (g = go) and exponential atmosphere (p = po e ' ^ y ), the equations of motion 
can be solved for velocity V and acceleration dV/dt (Allen and Eggers, 1958) 


V = V e 

E 


■ / (2pBPsiny ) 


dV / dt = E_ e “P y e 

2BP 


-Py / 

/ /(PBPsiny E ) 
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Runge-Kutta-Fehlberg routines are available in commercial software application packages designed for scientific and 
engineering computations. 
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where 


V E = entry velocity, m/s 

Ye = entry angle (down is negative) 

e = 2.7182818 (base of natural logarithms) 

P = inverse scale height (exponential density approximation) 
po = 1.2250 kg/m 3 
y = geometric altitude, m. 

Maximum deceleration - dV/dt| max occurs at altitude y = (1/p) In [-po/(P BP sin y E )] and is a function 
of entry velocity V E and flightpath angle at entry y E 

-dV /dt l MAX = pV E 2 siny E /2e 

At peak deceleration 

V = e~ 1/2 V E * 0.61 V E 

Maximum deceleration for “steep” nonlifting entry is presented as a function of flightpath angle at 
entry and entry velocity in table 10 and figure 20. Extremely high decelerations (greater than 
250 Earth g) are encountered for very steep entry angles (y E < -50 deg). 


TABLE 10. DIMENSIONLESS MAXIMUM DECELERATION (-dV/dt)/g FOR 

“STEEP” NONLIFTING ENTRY 3 


Entry 

Velocity, V E , 
km/s 

Flightpath Angle at Entry, y E , deg 

-10 

-30 

-50 

-70 

-90 

11 

54 

158 

241 

296 

315 

12 

65 

188 

287 

353 

375 

13 

76 

220 

337 

414 

440 

14 

89 

255 

391 

480 

511 

15 

102 

293 

449 

551 

586 


a g = 9.80665 m/s 2 and 1/p = 7200 m. 

The altitude at which peak deceleration occurs is a function of the ballistic parameter and flightpath 
angle at entry; it is independent of entry velocity. Peak deceleration altitude decreases as ballistic 
parameter increases and as the flightpath angle at entry decreases. 33 Note that slender bodies (lower 
drag) usually have higher ballistic parameters than blunt bodies. The altitude for maximum 

33 This reference is to the actual value of the flightpath angle at entry, not to its absolute value. (A change of entry angle 
from -30 deg to -45 deg represents a decrease in entry angle.) 
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deceleration for “steep” nonlifting entry is presented as a function of flightpath angle and ballistic 
parameter in table 1 1 and figure 2 1 . 


TABLE 11. ALTITUDE IN KILOMETERS AT MAXIMUM DECELERATION FOR 

“STEEP” NONLIFTING ENTRY" 


Ballistic 
Parameter, BP 
kg/m 2 

Flightpath Angle at Entry, y E , deg 

-10 

-30 

-50 

-70 

-90 

5 

66.4 

58.8 

55.7 

54.3 

53.8 

10 

61.4 

53.8 

50.7 

49.3 

48.8 

25 

54.8 

47.2 

44.1 

42.7 

42.2 

50 

49.8 

42.2 

39.2 

37.7 

37.2 

100 

44.9 

37.2 

34.2 

32.7 

32.3 

200 

39.9 

32.3 

29.3 

27.7 

27.3 


' p 0 = 1.2250 kg/m 3 and 1/p = 7200 m 


4.5 Lifting and Shallow Entry 

Chapman found a simple and more general solution to the equations of motion by introducing a 
coordinate transformation (Chapman, 1959). Numerical solutions have been tabulated for the 
following conditions (Chapman and Kapphahn, 1961): 

8 km/s < V E < 24 km/s 

-70 deg < y E < 0 

-4 < L/D < 4 


Chapman and Kapphahn’s solutions can be used for nonlifting entry (L/D = 0) and for lifting entry 
(negative as well as positive lift-drag ratios). 

Chapman was the first U.S. scientist to study the atmospheric entry corridor problem (see Chapman, 
1960). If a ballistic entry vehicle enters the atmosphere at too shallow an angle, it will skip out of the 
atmosphere. If the flightpath entry angle is too steep, the vehicle may burn up or suffer damage due 
to excessive structural loads. Guidance accuracy requirements are also discussed in Chapman’s work, 
Temperature constraints are discussed in section 6 (see Kolodziej, 1997). 
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4.6 Detailed Entry Trajectory Calculations 

Detailed trajectory simulations can be perfonned using three- or six-degree-of-freedom trajectory 
codes. The Program to Optimize Simulated Trajectories (POST) code has been used for this purpose 
for many years. POST does three degree-of-freedom (DOF) simulation. There is also a six-degree-of- 
freedom version called 6D POST. 34 POST has been used successfully to solve a wide variety of 
atmospheric entry problems. The program is capable of simulating and optimizing trajectories for 
space vehicles. 

Trajectory simulation codes include various atmospheric and gravity models. Aerodynamic 
coefficients (Newtonian, free molecular, flight- or ground-test data or CFD results) and aerodynamic 
heating models (both laminar and turbulent) are usually included as well. Some codes include 
propulsion models and reaction control system models. Navigation, guidance and flight controls can 
be simulated as well. 

Trajectory simulations are used at various stages of the design process. Generally, this is 
accomplished by selecting: a) appropriate atmosphere and gravity models; b) vehicle size, shape, 
mass, center of mass, TPS material, propulsion model; and c) initial conditions (entry coordinates, 
velocity, flightpath angle, and vehicle attitude). 

The following trajectory output results are obtained: a) altitude, velocity and acceleration profiles; b) 
vehicle attitude and stability; c) aerodynamic loads, M, Re, and dynamic pressure; d) heat transfer 
rates and heat loads; e) vehicle thermal response (heatshield and structure temperatures); and f) 
impact location. 35 

Results are evaluated to determine if mission requirements and constraints are satisfied. If not, 
vehicle size, shape, and material characteristics, as well as initial conditions, are modified until a 
satisfactory design is obtained. 

A trajectory simulation code called TRAJ (Allen et ah, 2005) can be used for 3- and 6-DOF 
trajectory simulation and thennal protection system design. 36 The TRAJ code has been validated with 
flight and ground test data for both trajectory variables and heatshield thennal response. Examples of 
TRAJ results for Stardust entry are presented in sections 6 and 8 (see figs. 30 through 34). 


34 At one time many software packages developed by government agencies and by government contractors were available 
to industry and to the public at cost from COSMIC or from NASA's Software Technology Transfer Center. Many of 
these items can be obtained from Open Channel Software (Open Channel Foundation, 2006). The POST code is 
described in NASA contractor reports (Brauer et al., 1975). 

35 Monte Carlo statistical simulations can be used to determine downrange and crossrange impact dispersions. This 
requires knowledge of variations in SRV properties, dimensional tolerances, accuracy of approximations, atmospheric 
temperature profdes, winds, and several other variables. It is necessary to run a very large number of trajectories to get 
accurate results. 

36 Copies of the TRAJ code are available from NASA with approval of an export control official. Point of Contact: Gary 
A. Allen, Jr., gallen@mail.arc.nasa.gov, 1-650-604-4228, MS 230-3, NASA Ames Research Center, Moffett Field, CA 
94035-1000. 
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4.7 Stability 


Vehicle stability during atmospheric entry is not a simple matter. The vehicle flow field is 
complicated and involves several flow regimes (from free molecular to subsonic). It is often 
necessary to perfonn full six-degree-of-freedom (6-DOF) trajectory simulations even for simple 
vehicle geometries. Six-DOF trajectory simulations require knowledge of stability derivatives, mass 
properties (weight and CG locations) and moments of inertia. Stability derivatives are usually 
determined experimentally. An additional equation of motion is required to account for rotation: 

19 = SM = qSLEC m 


where 

q = dynamic pressure 
S = reference area 

L = reference length (for pitching moment coefficient) 

I = pitch moment of inertia about the center of gravity 
9 = vehicle pitch angle measured from local horizontal 
M = pitching moment 

C m = pitching moment coefficient = M / (q S L). 
and Z indicates the sum of the moments about the center of gravity. 

In the hypersonic speed range the vehicle is usually statically stable (pitch oscillation amplitude 
decreasing as altitude decreases) if the static margin is adequate (center of gravity forward of center 
of pressure). However, it may be necessary to spin the vehicle to provide additional dynamic 
stability. The vehicle may also experience coning and tumbling motions. 

4.7.1 Static Stability 

For steady, horizontal flight, the pitching moment M about the vehicle center of mass (see fig. 12) is 

M = - N (x c r-x cG ) 

C m qSL = — C N qS(x cp -x CG ) 

Cm a = -CN a (X cp - X CG ) / L 
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where 


C N = N/qS 


C % = — — = pitching moment curve slope 
5a 

C N = = normal force curve slope. 

a 0a 

The vehicle is stable if C m<x is less than zero and if the center-of-gravity location xcg is forward of the 
center-of-pressure location x C p: 


C m<x < 0 and x CP - x CG > 0 

The quantity (x G p - x G g) is often referred to as the static margin. The larger the value of the static 
margin the more stable the vehicle in the sense of static stability (Etkin, 1972). 

A sharp, conical body in steady horizontal, flight will be statically stable in the pitch plane if C m(X is 
negative and unstable if C ma is positive. The Newtonian approximation for the pitching moment 
curve slope for sharp, conical bodies is shown in figure 22. Figure 22 also shows the effect of CG 
location x G g- A sharp conical body becomes less stable as the CG location moves aft (closer to the 
downstream-facing base of the cone). For a sharp cone at zero angle of attack, the CP location x G p is 
given by Newtonian theory 


x C p = 2/3 F ( 1 + tan 2 9c) 

where F = cone length and 0 G = cone semivertex angle. 

4.7.2 Dynamic Stability 

In atmospheric entry the motion of an entry vehicle may not be stable. As altitude decreases dynamic 
pressure first increases to a maximum and then decreases to impact (or until terminal velocity is 
reached). The increase in dynamic pressure results in pitch oscillations if the vehicle axis is not 
aligned with the velocity vector at entry. 

Tobak, Allen, and Sommer studied hypersonic dynamic stability for atmospheric entry at Ames in the 
1950s. They established dynamic stability criteria and analytically determined the amplitude and 
frequency of pitch oscillations (Tobak and Allen, 1958 and Sommer and Tobak, 1959). They also 
derived an expression for pitch angle as a function of altitude in terms of Bessel functions 
(Allen, 1957). 

An application of Tobak’s method is to determine the dynamic stability parameter as a function of 
altitude and estimate the convergence-divergence boundary for unstable pitch oscillations (i.e., the 
altitude at which the sign of the dynamic stability parameter changes from positive to negative). 

Table 12 shows the results of this calculation for three flightpath entry angles. In each case pitch 
oscillations due to finite angle of attack at entry converge at altitudes above the convergence- 
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divergence altitude shown in table 12 and diverge at lower altitudes. Note that these results are valid 
only for “steep” entry angles and hypersonic speeds. 


TABLE 12. CONVERGENCE-DIVERGENCE ALTITUDE FOR PITCH OSCILLATION 
STABILITY FOR “STEEP” NONLIFTING ENTRY (BP = 50 kg/m 2 , V E = 15 km/s AND 

y E = 125 km) a 


Flightpath angle at entry, y E , deg 

Convergence-divergence altitude, km 

-20 

44.9 

-40 

40.4 

-60 

38.3 


a For L/D = 0 and M » 1. Note that these results cannot be extrapolated to altitudes below 34 km. 
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5. DECELERATION AND RECOVERY 


In general, deceleration and recovery planning begins with trajectory computations (see sec. 4). The 
simplest case is nonlifting entry where the critical variable, the ballistic parameter 37 BP is constant. 
More complicated cases arise when constant or variable lift is introduced to alter the trajectory. Lift 
can be used to change orbit plane, reduce heat flux or total heat load, or improve targeting accuracy. 
It is also possible to control the aerodynamic drag force acting on the vehicle; this is called “drag 
modulation.” 


5.1 Drag Modulation 

Drag modulation can be used to change the speed of a nonlifting vehicle. Drag increases as the 
vehicle cross-sectional area (in a plane normal to the velocity vector) increases. If lift is assumed 
equal to zero, increasing or decreasing the vehicle cross-sectional area can modulate drag. Drag 
modulation can be used to constrain the trajectory variables and control the impact location. 
However, drag modulation is less effective than varying the lift of the vehicle (lift modulation). 

A drag modulation concept was proposed for an Aero- Assist Orbital Transfer Vehicle (AOTV) in the 
1980s. This concept was a combination of a balloon and a parachute and termed a “ballute.” Unlike 
wake-deployed ballutes the AOTV ballute inflated around the vehicle. Drag modulation would occur 
by increasing or decreasing the pressure inside the balloon, which would in turn change the cross- 
sectional area of the vehicle. The concept attracted much interest due to its apparent simplicity. 
However, it was difficult to implement and the concept was eventually dropped. The major difficulty 
was providing thermal protection for the inflatable structure. The Soviet Union has used inflatable 
braking devices for planetary entry in the past and Russia continues to study their use. 


5.2 Parachute Recovery 

For entry vehicles with low values of L/D at subsonic speeds the use of a parachute for final descent 
and landing is an attractive option. This approach has been used for blunt entry capsules (both Earth 
and planetary; Apollo and Mars Pathfinder, for example) and also proposed for lifting vehicles with 
poor subsonic flight qualities. Mach number and dynamic pressure, along with Reynolds number, are 
the critical variables. Lower Mach numbers are preferred as it is difficult to deploy a parachute at 
high supersonic Mach number and high dynamic pressure. The problems include: deployment shock 
to the material; uncertainty of proper parachute inflation; and difficulty of predicting overall 
performance (particularly in the transonic-supersonic range where shock waves and uncertain 
pressure distributions are encountered). 


37 The ballistic parameter BP is defined here as m/Cn S and its units are kg/m 2 . A similar quantity in English units is 
called the ballistic coefficient and is defined as W/ C D A and given in units of lbm/ft 2 , where W is the body weight and 
the reference area A is equivalent to the area S in the definition of ballistic parameter. 
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For most applications, parachute deployment usually occurs around Mach 1.8 because of dynamic 
instability at lower Mach numbers. Since the flow field is greatly influenced by body shape ground or 
flight testing is usually required. 


5.3 Final Descent and Landing 

A heatshield having significant thermal mass can transfer thermal energy to the interior of a SRV. 
Even though the heatshield surface temperature decreases rapidly as velocity decreases at lower 
altitudes heat conduction from the heatshield to the vehicle interior can continue after landing. For 
Space Shuttle Orbiter flights structural temperatures beneath the TPS tiles would continue to increase 
after the vehicle had landed if external cooling was not applied. 

During the final descent of the Mars Pathfinder the forward heatshield was separated from the lander 
vehicle. This method could also be employed on a Mars SRV entering the Earth’s atmosphere. The 
Mars Pathfinder deployed a drogue chute and air bags that inflated just prior to impact. 


46 



6. HEAT TRANSFER AND THERMAL PROTECTION MATERIALS 


6.1 Heat Transfer 

Heat transfer prediction is critical to the design of a thennal protection system (TPS). Heat flux, total 
heat load and surface temperature histories affect the selection of: an appropriate heat-shield 
material; the heat-shield thickness; and the means of attaching the heatshield to the vehicle structure. 
In the 1950s, it was thought that there was a “thermal barrier,” much like the belief that there was a 
“sonic barrier,” that was impossible to go beyond. It was believed that the very high kinetic energy 
associated with an entry vehicle would greatly exceed the heat of vaporization of most materials. 
However, only a small fraction of the vehicle kinetic energy goes into actually heating the vehicle. 
Most of the energy is used to heat the air ahead of the vehicle as it is compressed by the bow shock 
wave. The thennal barrier may be readily penetrated with appropriate TPS design. (Allen’s discovery 
that bluntness reduces the effect of entry heating is discussed in sec. 4.2.) 


6.2 Radiative Heat Transfer 

The dominant form of entry vehicle heat transfer is convective heat transfer at the vehicle surface. In 
some situations radiative heating from the shock layer gasses can be very significant. At vehicle 
velocities above 13 km/s in the earth’s atmosphere, the shock layer becomes extremely hot (see 
fig. 2). At these velocities radiation from the shock layer to the vehicle surface can become a 
significant fraction of the total incident energy. Computational methods for radiative heat transfer 
effects are extremely difficult and become more difficult with the introduction of ablation products 
into the boundary layer. An added complexity is that thennodynamic and transport properties are not 
well understood at temperatures above 9000 K. 

Flight testing may be necessary to detennine radiative heat flux at conditions that would be 
experienced by an SRV. In some situations simple approximations for radiant heat flux qRAD can be 
used. At the stagnation point on the body: 

q rad t = C r N A p B f(V) 

where: A is a function of p and V; B and C are constants; and f(V) is a function of velocity for a 
given atmosphere (Tauber and Sutton, 1991). 

Radiative heating was a significant portion of the total heat flux for the Galileo Jupiter entry probe. 
For the Apollo Command Module approximately 34 percent of the total heat flux was due to 
radiative heating (see table 13). For a typical nonlifting entry body as much as 83 percent of the total 
heat load can be due to radiative heating (see table 13). 
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TABLE 13. PEAK STAGNATION POINT HEAT FLUX AND TOTAL HEAT LOAD FOR 

EARTH ENTRY SPACECRAFT 


Entry vehicle 

BP, 

kg/m 2 

V E , 

km/s 

Ye, 

deg 

qMAX t, 

W/crn 2 

Q, 

J/crn 2 

Apollo CM, L/D « 0.3 

« 500 

* 11 

— 

« 5 1 0 b 

— 

a 

Stardust 

68.2 

12.9 

-8.2 

856 

23,730° 

Entry Body, L/D = 0 

50 

15 

-10 

1750 

19,340 d 

Entry Body, L/D = 0 

50 

15 

-50 

11,170 

20,300° 


a Results of NASA TRAJ simulation, 1999. 

b Approximately 34 percent of the total heat flux is radiative and 66 percent is convective. 
c The total heat load is 9 percent radiative and 91 percent convective. 
d The total heat load is 61 percent radiative and 39 percent convective. 
e The total heat load is 83 percent radiative and 17 percent convective. 

The TRAJ simulation results for a body shaped like the Stardust entry capsule is shown in figures 23 
and 24. At entry velocities less than 12 km/s the stagnation point heat flux and the stagnation point 
heat load are dominated by convective heating. As the entry velocity increases the radiant heat flux 
increases rapidly. At Ve = 15 km/s, the peak stagnation point radiative heat flux is more than two 
times greater than the peak stagnation point convective heat flux. At entry velocities greater than 
15 km/s, the heat load due to radiation is greater than the heat load due to convection. 


6.3 Heat Transfer at Vehicle Surface 

At the surface of an entry vehicle heat is transferred from the hot gases in the shock layer to the 
vehicle surface by forced convection qcoNv and through the vehicle structure to the vehicle interior 
by conduction qcoND- The convective heat flux can be estimated using the methods given in 
section 3.8. Methods for estimating conductive heat flux are given in textbooks on heat transfer 
(Jacob, 1949). The Stefan-Boltzmann radiation law gives the re-radiative heat flux 38 at the vehicle 
surface: 


where 


qRE-RAD 


S s a (Tw 4 - Tspace 4 ) 


S = surface area 

s = emissivity of material surface 

a = Stefan-Boltzmann constant = 5.67040 x 10 8 W/(nr-K 4 ) 
T w = wall temperature 

Tspace = temperature of space or radiation sink. 


38 This is the heat transferred from the vehicle surface to the space environment by radiation. 
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Note that the fourth power of T w shows that re-radiation is far more significant at high temperature. 
For practical considerations, it is desirable to minimize heat conduction to the vehicle interior (which 
is why low-density ceramics are often used) and to maximize the heat re -radiation (which is why 
surface coatings with high emissivity are desirable). Values of emissivity for a variety of TPS 
materials are given in NASA’s TPSX online database (Squire, 2006). 39 


6.4 Heatshield Design Considerations 

In designing a heatshield for an atmospheric entry vehicle, there are a variety of factors to consider. 
Two dominant factors are system weight and backwall design temperature; the latter determines how 
much thermal energy will flow to the vehicle interior. For large reusable vehicles, the economics of 
the system, as well as qualification costs, are of considerable concern. As an example, the Chinese 
have successfully used oak wood for an entry capsule TPS. The oak performed well and is 
economical but results in a fairly heavy vehicle. Lightweight ceramic materials have much lower 
thennal conductivity and much lower density than oak, but are much higher in cost. The payoff for 
lightweight TPS materials is the added mass that can be used for scientific instruments or other 
payload. SRV missions have large propellant mass requirements. If TPS mass can be reduced, less 
propellant will be required (or payload mass can be increased). 

A method for heat shield design can be described as follows. First, the mission requirements are 
defined and trajectory calculations are made for specific entry conditions. Then first order 
calculations are made for critical heating areas on the vehicle at points along the flightpath. Thermal 
analyses are perfonned to determine required TPS thicknesses at critical locations. For TPS materials 
attached to the vehicle structure with high-temperature adhesives, the backwall temperature cannot 
exceed the maximum allowable temperature of the adhesive. Different materials can be selected for 
different locations on the vehicle. For example, rigid ceramic tiles can be used on windward (or 
leading edge) surfaces and flexible blanket materials such as advanced flexible reusable surface 
insulation (AFRSI) or tailorable advanced blanket insulation (TABI) can be used on leeward 
surfaces. At a later stage, a more detailed flowfield analysis of the vehicle can be performed to 
achieve optimal sizing. Analytical and numerical methods have been developed by NASA to greatly 
expedite this design process. Various trajectory simulation codes can be used for TPS sizing. 


6.5 Methods of Thermal Protection for Sample Return Vehicles 

Thermal protection system design is critical for sample return vehicle mission success. SRVs 
entering the Earth’s atmosphere at velocities greater than 10 km/s will encounter severe heat fluxes 
and heat loads as shown in table 13. Methods for estimating convective and radiative heat fluxes for 
atmospheric entry vehicles are discussed in section 3. 

An effective thennal protection system is required to assure survival of sample return capsules and to 
ensure that capsule interior temperature limits are not exceeded. An insulating material may be used 
to protect the load-bearing structure or an active cooling system may be employed. In severe heating 


39 


Registration is required; some restrictions apply. 
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environments the insulating material may be allowed to ablate. Some common thennal protection 
schemes are described below. 

6.5.1 Heat Sink 

A heat-sink structure is a simple, passive thermal protection method. 40 It can be a homogeneous or 
composite material or a honeycomb structure. Composites are made up of two or more constituents 
bonded together to optimize thermal and/or mechanical properties. Composite materials are generally 
anisotropic (properties have different values when measured in different directions). In many 
applications heat-sink materials provide inadequate strength, excessive weight or uncertain 
survivability when subjected to entry heating. (Beryllium-copper was used as a heat-sink material for 
early re-entry vehicles, but fabrication was very difficult.) 

6.5.2 Radiative Surfaces 

Radiative surfaces dissipate a portion of the incident radiant energy by re -radiation. If the structure 
can tolerate radiation equilibrium temperatures, which vary as q 14 , and the material has a low thermal 
conductivity, then radiative cooling can be very effective. The ratio of emitted radiant flux density to 
incident radiant flux density is called emissivity s. 

6.5.3 Film and Transpiration Cooling 

Film or transpiration cooling can be used to protect vehicles from high entry-heating environments. 

In film cooling, a fluid is injected into the boundary layer near the stagnation point to protect the 
surface. In transpiration cooling, the fluid is injected through pores in the material surface. Film and 
transpiration cooling are less effective when the heat input is predominantly radiative rather than 
convective. Analyses and experiments show that the heat flux can be reduced significantly by 
increasing the mass addition rate (Tauber, 1989). The effectiveness of film or transpiration cooling 
can be further increased by selecting a coolant gas with a molecular weight lower than that of the 
boundary layer gases. Film and transpiration cooling are more effective in laminar flow than in 
turbulent flow. A disadvantage of film and transpiration cooling is that a complex supply and 
regulating system may be required. However, film or transpiration cooling may be required to protect 
small areas with sensors or windows. 

6.5.4 Ablation 

Ablation cooling is a very common method of thennal protection for planetary entry vehicles 
(e.g., Apollo Command Module, Mars Pathfinder and Galileo). Ablation materials are consumed as 
they are heated to sublimation temperature. Sublimation phase change (solid state to liquid or vapor) 
and gaseous transpiration promote surface cooling (material surface temperature affects the amount 
of heat radiated from the surface). Ablation is a complex process and can involve a variety of 
chemical reactions including: combustion; sublimation; vaporization (liquid to gas); erosion; 
oxidation; charring; pyrolysis; melting; molecular dissociation; and recombination. Figure 25 is a 
schematic diagram of the charring ablation process. 


40 A heat sink is a body that can absorb or reject a large quantity of heat without an appreciable change in temperature. A 
heat sink is also referred to as a heat reservoir. 
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NASA defines ablation as (Allen, 1965): The removal of surface material from a body by 
vaporization, melting, chipping, or other erosive process; specifically the intentional removal of 
material from a nose cone or spacecraft during high speed movement through a planetary 
atmosphere to provide thermal protection to the underlying structure. 

The following description appears in the current online version of the NASA’s Dictionary of 
Technical Terms (Glover, 2006): Ablating materials are used on the surfaces of some reentry 
vehicles to absorb heat by removal of mass, thus blocking the transfer of heat to the rest of the 
vehicle and maintaining temperatures within design limits. Ablating materials absorb heat by 
increasing in temperature and changing in chemical or physical state. The heat is carried away from 
the surface by a loss of mass (liquid or vapor). The departing mass also blocks part of the convective 
heat transfer to the remaining material in the same manner as transpiration cooling. 

As shown in figure 25, vaporization and sublimation occur at the outer surface of the ablation 
material. Oxidation and combustion also occur at the outer surface and result in surface recession as 
the ablation material is consumed. The ablating surface may serve as a catalyst to promote 
recombination of oxygen and nitrogen atoms in the boundary layer. Pyrolysis, a chemical change 
caused by heat addition, begins as thennal energy penetrates into the material (at around 670 K for 
phenolic resins). About one-half of the material volatizes and the other half is left as a carbonaceous 
residue or char. Decomposition gases from the pyrolysis zone percolate through the char layer and 
absorb thennal energy. 

Wood is an example of a charring ablation material. Synthetic resins cured with reinforcing fibers are 
used as charring ablators for entry vehicle heat shields. The energy absorbed per unit mass ablated is 
a measure of ablation effectiveness. For silicon nitride, the energy absorbed per unit mass ablated is 
6 MJ/kg (Tauber, 1989). Properties of typical TPS materials are given in table 14. 

The carbon-phenolic materials used on the Galileo probe may be too heavy (1200 to 1920 kg/m ) for 
SRVs. The heat shield for the Apollo command module was made from AVCOAT-5026 (529 
kg/m ). Superlight ablators (SLAs) were used on the Mars Viking spacecraft in 1976 and the Mars 
Pathfinder spacecraft in 1997. SLA-561V has a density of 265 kg/m\ A disadvantage of SLA 
materials is their high manufacturing cost and fragility. Phenolic impregnated carbon ablator (PICA) 
is a lightweight ceramic ablator (LCA) material. PICA can be made in densities from 225 to 880 
kg/m 3 and was selected for the Stardust sample return mission. 41 

Ablation material sizing can be perfonned using the Fully Implicit Ablation and Thermal Response 
(FIAT) code developed by Y.-K. Chen (Chen and Milos, 1997). The FIAT code is incorporated in 
NASA’s TRAJ trajectory simulation code (Allen et al., 2005). The code is easy to use, but only four 
ablation materials are available. The user first specifies the vehicle geometry and initial conditions 
(y E , Ve, and Ye) and selects the heatshield material and thickness, then runs the trajectory simulation. 
The code calculates heatshield mass, surface and in-depth temperature-time histories, surface 
recession and consumable mass flux for the calculated trajectory. Another option is to allow the code 
to calculate an optimal heatshield thickness for a specified backface temperature. The second option 
does not always provide solutions. Typical results for a Stardust sample return mission are given in 


41 PICA is composed of a fibrous carbon substrate (manufactured by Fiber Materials, Inc.) and a phenolic transplant. 
PICA shows high effective heat of ablation at heat fluxes between 300-1500 W/cm 2 . 
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table 15. Typical results for a nonlifting sample return entry vehicle are given in table 16. (Note that 
these results were obtained using a 1999 version of the TRAJ code.) 


TABLE 14. PROPERTIES OF SOME TPS MATERIALS* 


Material 

Density, 

kg/m 3 

Thermal 

Conductivity, 3 

W/m-K 

Emissivity, 3 

8 

Melt 

Temp., 

K 

A VC OAT 5026 (virgin) 

529 

0.242 

0.67 

1920 

A VC OAT 5026 (char) 

264 

— 

0.49 

— 

Carbon-Phenolic 
(Acusil I) 

480 

0.112 

0.95 

2000 

Reinforced Carbon- 
Carbon (RCC) 

1580 

5.05 (normal) 
7.88 (parallel) 

0.54 to 0.9 

2030 b 

Teflon 

2190 

0.251 

— 

— 

Shuttle Tile 
(LI-900) 

144 

0.047 

0.88 

1755 b 

Beryllium 

1840 

170 

— 

8 1 0 b 

Oak 

610 

0.146 

— 

— 


* Squire, 2006 (TPSX database). 

a Thermal conductivity and emissivity are functions of temperature. 
b Single-use temperature limit. 


TABLE 15. TRAJ PLUS FIAT ABLATION SIMULATIONS FOR STARDUST SRV 


Material 

Q, 

kJ/ctn 2 

Qrad, 

kJ/ctn 2 

Heatshield 

Thickness, 

cm 

Heatshield 

Mass, 

kg 

Surface 

Recession, 

cm 

Mass 

Loss, 

kg 

PICA 3 

27 

2 

5.08 

— 

1.07 

1.06 

PICA-15 b 

23.7 

2.1 

5.08 

7.55 

0.76 

0.26 c 

Carbon- 

Phenolic b 

23.7 

2.1 

3.81 

33.97 

0.12 

— 


11 FIAT simulation (Chen and Milos, 1998). Bondline (ablation material backface) temperature < 523 K. 
b TRAJ plus FIAT simulation, 1999. Initial conditions: y E = 132.9 km, V E = 12.86 km/s, y E = -8.2 deg for 
BP = 68.2 kg/m 2 . 

c Mass loss estimated from integral of char plus gas mass flux. 
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TABLE 16. TRAJ PLUS FIAT ABLATION SIMULATION FOR A NONLIFTING SRV 


Material 

Q, 

kJ/cm 2 

Qrad, 

kJ/cm 2 

Heatshield 

Thickness, 

cm 

Peak 
Surface 
Temp., K 

Bondline 
Temp., K 
at t = 94 s 

Surface 

Recession, 

cm 

PICA- 15 

28.0 

4.4 

3.43 

3536 

440 

1.36 

Carbon- 

28.0 

4.4 

2.54 

3262 

439 

0.18 

Phenolic 








a Initial conditions: y = 125 km, V E = 15 km/s, y E = -10 deg for BP = 50 kg/m 2 (Blunt cone: 9 C = 25 deg, 
r N = 0.171 m, r B = 0.43 m). Initial bondline (ablation material backface) temperature = 239 K. 


6.6 Ablation Constraint Determination 

Kolodziej has developed a spreadsheet method to calculate atmospheric-entry constraints (Kolodziej, 
1997). If vehicle characteristics are specified, a bounding altitude-velocity curve can be obtained that 
indicates a limit for a system variable (such as stagnation point temperature). For example, if it is 
necessary that the stagnation point temperature not exceed the material ablation temperature, then a 
bounding curve can be determined that separates the region where ablation occurs from the no- 
ablation region as shown in figure 26. For this example ablation occurs when the trajectory curve 
falls below or to the right of the “ablation-constraint” curve. This method includes hot wall, wall 
catalysis, imperfect gas properties and rarefied flow effects. An example is presented below. 

For materials that are good insulators and radiate efficiently: 

_ rj, 4 

qcoNv - qcoND - a s t w 

If qcoNv » qcoND, h w « h t and T w = T M ax (the maximum nonablating use temperature of the 
material), then for an axially-symmetric body with nose radius tn 

_ t 4 

qcoNv - a s t max 


and 

qcoNv = const p 17 " V 3 / rN 17 “ 


or solving for V = Vconstraint 

V constraint = (o s/0.000183) 173 Ov/p) 176 Tmax 4/3 , m/s 

If freestream density p = po e' Py (where p 0 = 1.225 kg/nr and 1/p = 7200 m), then Vconstraint is a 
function of altitude alone for given values of s, T M ax, and r N . Ablation-constraint curves are shown in 
figures 27 through 29. The curves show the effects of emissivity, ablation temperature and nose 
radius on the ablation constraint. These curves correspond to values that might be expected for a 
noncatalytic cold wall. For a noncatalytic hot wall, the method shows the constraint is less 
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conservative (the curve is shifted down and to the right). For fully catalytic walls the constraint curve 
falls between the hot-wall noncatalytic and the cold-wall noncatalytic curves. 


6.7 Hypersonic Simulation and Optimization Codes 

The Hypersonic Vehicle Optimization Code (HAVOC) can be used to obtain estimates of TPS 
requirements (thickness, weight and cost) for a specified TPS material, vehicle geometry and 
trajectory. The HAVOC engineering design code has been used for sizing of nonablative TPS 
materials for hypersonic aircraft (Kolodziej et ah, 1998). 

An aerodynamic heating and TPS sizing code (MINIVER) can be used both to predict the 
aerothermal environment and to perform simple TPS sizing for aerospace vehicles that operate in the 
hypersonic flight regime (Wurster et al., 1999). The code uses approximate heating methods for 
predicting heat flux and perfect-gas or equilibrium-air chemistry for aerodynamic flow field 
approximations. 

MINIVER has been used as design tool in government and industry for a wide variety of vehicle 
configurations. The code has the ability to quickly provide the thermal environments required for 
TPS sizing. MINIVER was originally developed at McDonnell Douglas in the 1970s under 
government contract. Many government, military, educational, and aerospace industry installations 
currently use it. 42 

6.7.1 Aerothermodynamics 

Analytical methods are used to detennine aerodynamic coefficients for various geometric 
configurations. Experimental data can be used to improve the results. In the supersonic/hypersonic 
speed range, lift and pressure drag are computed using tangent wedge/tangent cone methods. 
Hypersonic pressure drag is found using Newtonian theory. HAVOC includes base pressure drag, 
skin friction drag, and approximations for heat transfer coefficients. At present MINIVER does not 
include provisions for radiative heat transfer. 

6.7.2 Structural Analysis 

The HAVOC code has structural analysis and weight estimation capabilities. After body stress 
distributions are calculated, the minimum amount of structural material required to prevent failure is 
determined. Tensile yield, compressive yield, local buckling, and gross buckling failure modes are 
considered. TPS thickness determination is based on quasi-steady state, one-dimensional heat 
conduction analysis and the requirement to keep the interior of the vehicle below a specified 
temperature. 


42 NASA Point of Contact: Kathryn Wurster, NASA Langley Research Center, Mail Stop 451, Hampton, VA 23681- 
2199, 1-757-864-4487, K.E.Wurster@nasa.gov. 
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6.7.3 Trajectory Simulation 

A three-degree-of-freedom trajectory simulation is available in the HAVOC code or a Mach number 
versus altitude profile can be specified. The equations of motion are for a point-mass body moving 
relative to a rotating, spherical earth. The trajectory can be constrained by maximum dynamic 
pressure or by maximum structural temperature. 
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7. GROUND AND FLIGHT TEST SIMULATION 


A variety of analytical and experimental methods are available for use in SRV design. These methods 
may include a combination of ground tests to simulate portions of the flight regime or aspects of flow 
phenomena (e.g., aerodynamic loads, TPS thermal response or chemical reactions). As the flight 
concept matures simplified flight experiments can be perfonned to add important design information 
that is not available from ground facilities or computation (e.g., boundary layer transition at 
hypersonic velocities). The final validation is of course the actual flight of the SRV. Some of the 
critical simulation parameters (Mach number, Reynolds number, Knudsen number, etc.) were 
discussed in section 1 . CFD methods are used extensively to detennine aerodynamic and 
aerothennodynamic flow characteristics. However, ground- and flight-testing are still necessary to 
validate and “calibrate” the CFD methods and their results. 43 

An inventory of major government aerospace test facilities is available on the Internet (NASA 
Headquarters, 2006). 


7.1 Ground Test 

Aerodynamic ground testing facilities and techniques are described by Lukasiewicz (1973) and by 
Pope and Goin (1965). Brief descriptions of the principal types of ground-test facilities are given 
below. 

7.1.1 Conventional Wind Tunnel 

Conventional wind tunnels are used to determine aerodynamic forces and moments as well as vehicle 
surface pressure distributions. Compressibility and viscous effects are simulated by matching test 
Mach number and Reynolds number with flight values. For very high altitudes, Knudsen number is 
matched in low-density wind tunnels. Strain gauge balances are used to measure forces and moments. 
Pressure transducers are connected to pressure ports to measure surface static pressures. Nonintrusive 
velocity measurements are made with a laser Doppler velocimeter. Hot-wire anemometers are used to 
measure turbulent velocity fluctuations. There are a variety of flow visualization techniques to 
determine shock wave shape and location, surface streamline patterns and heat transfer 
characteristics 44 (Liepmann and Roshko, 1957). 

At high Mach numbers (above around Mach 7) conventional wind tunnels are limited by test section 
size and run time. It is usually necessary to use small models and rapid-response data acquisition 
systems. There are also limitations on maximum available stagnation temperatures and pressures. 


43 In cases where there is little or no previous CFD experience, it is often necessary to vary the initial conditions and/or 
the turbulence model to find acceptable solutions; i.e., a solution that converges and matches reality. 

44 Schlieren, shadowgraph, and interferometer techniques have long been used in conventional wind tunnels for flow 
visualization of aerodynamic effects; e.g., shock wave and boundary layer effects. These methods are based on the 
optical property of a fluid that the speed of light varies with the density of the fluid through which it is passing. 
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7.1.2 Arc-Heated Facilities 


Arc heated facilities, or arc -jet wind tunnels, heat air (or another gas) to very high temperatures to 
simulate atmospheric entry heating. These devices employ high-voltage electrical connectors around 
a cylindrical tube. As cold gas moves through the column it is heated by an electrical discharge and is 
expanded through a nozzle into a test chamber. The test chamber may be continuously evacuated by 
a vacuum system in order to maintain the desired test conditions. 

Although arc -jet facilities have proven extremely valuable in terms of TPS materials selection and 
development, there are certain limitations. First, the chemical composition of the flow in flight is 
significantly different from that in an arc -jet facility. In flight, the gas is not dissociated or ionized 
upstream of the bow shock. Second, the upstream flow values are not constant; the flow approaching 
the test article is not necessarily uniform (flow velocity vectors may be divergent rather than 
parallel). Third, the Mach number and the Reynolds number are usually much lower than in flight 
making it difficult to simulate compressibility and viscous effects. Thus, the arc jet is an important 
tool for simulating high-enthalpy flows, but arc -jet facilities are poor wind tunnels in the classical 
sense. Performance values for NASA Ames Research Center’s arc-jet test facilities are given in 
table 17. 


TABLE 17. NASA AMES RESEARCH CENTER ARC-JET TEST FACILITIES 


Facility 

Test Gas 

Mach 

Number 

Test Article 
Size, cm 

Bulk 

Enthalpy, 

MJ/kg 

Surface 

Pressure, 

atm 

Convective 
Heat Flux, 
W/cm 2 

Aerodynamic Heating 

Air, 

4-12 

20 dia 

12-33 

0.001-0.125 

0.06-255 

Facility (AHF), 
20 MW 

Nitrogen 


66x66 




Turbulent Flow Duct, 

Air, 

3.5 

20x51 

3.5 -9.3 

0.02-0.15 

2-68 

12 MW 

Nitrogen 






Panel Test Facility 
(PTF) 20 MW 

Air 

5.5 

35 x35 

4.6-33 

0.0005-0.05 

0.6-85 

Interaction Heating 

Air 

5.5 -7.5 

61 x 61 

7-46 

0.0001 - 1.2 

0.6 - 749 a 

Facility (IHF), 75 MW 



46 dia 





a The Interaction Heating Facility can produce a radiative heat flux of 23 W/cm 2 . 


7.1.3 Shock Tunnel 

A simple shock tube consists of two chambers containing gases at different pressures separated by a 
diaphragm. When the diaphragm is ruptured the high-pressure gas flows into the low-pressure 
chamber. If the pressure difference across the diaphragm is large enough, a supersonic flow and a 
travelling shock wave are created. Shock tubes provide valuable experimental data for high-speed gas 
flows but are limited by chamber pressures (or pressure ratio) and by chamber length. The 
consequences are very short run times and an inability to achieve hypersonic Mach numbers. 
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A shock tunnel overcomes some of the limitations of the shock tube by attaching a second diaphragm 
and a hypersonic nozzle at the end of the low-pressure chamber. When the second diaphragm is 
ruptured, gas flows through the nozzle and into a test section and then into a vacuum reservoir. It is 
possible to simulate dissociation, ionization and free molecular flow effects in a shock tunnel. 
Perfonnance values for typical shock tunnels are: 


Velocity 

Reynolds number 
Mach number 
Simulated Altitude 
Run time 


2 to 5 km/s 

3.3 x 10 4 to 3.3 x 10 s per m 
6 to 24 
10 to 60 km 
4 to 20 millisec 


Note that run times are very short and flow velocities are much lower than Earth-entry velocities for 
Mars sample return missions. These values are for the U.S. Defense Department’s Calspan shock- 
tunnel facility in Buffalo, NY (NASA Headquarters, 2006). 

7.1.4 Ballistic Range Facilities 

Ballistic range facilities are used to obtain static and dynamic aerodynamic data. In a ballistic range a 
model is fired from a gun into a test gas (usually a quiescent gas). Test data are measured and 
recorded at stations along the model flightpath. Test Reynolds numbers can be varied by increasing 
or decreasing the ambient pressure of the test gas. Ballistic range testing is discussed in an AGARD 
document by Canning et al. (1970) and in an AIAA paper by Strawa et al. (1988). 

NASA Ames Research Center’s Hypervelocity Free Flight Facility has been used extensively for 
testing a variety of entry vehicle configurations. The facility is capable of launching a 22 mm dia. 
model at a speed of 8 km/s. See Chapman and Yates (1998) for a description of the use of a ballistic 
range in the design of a planetary probe. 


7.2 Flight Test 

7.2.1 General Considerations 

Ultimately in the design of an entry system, particularly one utilizing new materials or a new 
configuration, there is little substitute for an actual flight test. 45 This is the case if there are aspects of 
the flow field or configuration that are outside of the experimental database. Depending on the design 
aspect of specific interest (and of course the cost) different launch vehicles, flight velocities, or 
trajectories may be selected. 


45 The X-15 experimental aircraft achieved a maximum Mach number of 6.7 and a maximum altitude of 103.6 km in the 
1960s. Three rocket-propelled, air-launched X-15s were built. They were launched from a B-52 aircraft and flew 199 
flights from 1959 to 1969. 
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Note that in the National Aerospace Plane (NASP) program of the 1980s there was sufficient 
confidence in the CFD codes that experimental validation was considered unnecessary. However, it 
was later realized that flight testing was absolutely essential to “calibrate” and validate the 
computational models for boundary layer transition, engine inlet pressure distribution, etc. By that 
time, it was difficult to incorporate flight testing; this may have contributed to the subsequent 
cancellation of the program. 

A simple sounding rocket can reach Mach 8 (« 2.6 km/s) with relative ease. This allows meaningful 
aerodynamic experiments to be carried out at hypersonic speeds. For higher entry velocities, where 
aerodynamic heating becomes more significant, larger and more costly ballistic missiles may be used 
to achieve entry velocities on the order of 5 km/s. Most ICBMs, and many U.S. launch vehicles, fall 
into this category and have been extensively flight tested by the U.S. military. Boundary layer 
transition data were obtained on ICBM re-entry vehicle flight flights in the 1960s. High entry 
velocities can also be achieved using the “pile-driver” concept. The lower stages of the vehicle are 
used to gain altitude; then the final stage is pointed downward and its rocket engine ignited. The 
“pile-driver” method was used successfully in the NASA Ames Planetary Atmosphere Entry 
Technology (PAET) experiment — a four-stage Scout rocket launched from Wallops Island achieved 
an entry velocity of 7 km/s. 

For testing above 7 km/s (roughly equivalent to re-entry from Earth orbit), the vehicle may require 
even greater “pile-driver” thrust. The NASA-proposed Aeroassist Flight Experiment (AFE) of the 
late 1980s was to be deployed from the Shuttle Orb iter and use a solid-rocket motor to thrust the 
vehicle into the atmosphere at approximately 12 km/s. Other experiments in this velocity range were 
done for U.S. and proposed Soviet lunar return capsules (Apollo and Soyuz). The cost of a flight 
experiment usually scales directly with entry velocity. For typical flight test entry velocities see 
table 18. 


TABLE 18. TYPICAL FLIGHT TEST ENTRY VELOCITIES 


Vehicle/System 

Example 

Velocity, Ve, km/s 

Sounding Rocket 

Terrier/Black Brant 3 

2.5 

Sounding Rocket, Military ICBM 

Scout, Minuteman 

5 

Shuttle Orbiter and “pile driver” 

STS Orbiter, PAET 

7+ 

Blunt ablator and “pile driver” 

Apollo CM (unmanned) 

9 

Combination “pile driver” 

Project Fire b 

11.6 


a See Jane's Space Directory (Baker, 2005). 

b Project Fire was a NASA flight experiment with a blunt entry vehicle that achieved a velocity of 1 1 .6 km/s during 
re-entry (Dingeldein, 1965). An Atlas launch vehicle was used with and Antares rocket motor fired at the end of the 
coast phase. 
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7.2.2 Impact and Recovery 


Impact and recovery must be considered for every SRV experiment. Recovery phase design may 
have a significant bearing on how data is collected and analyzed and how accurately a probe or entry 
vehicle reaches its intended target. Significant issues affecting recovery are whether or not a 
parachute is used, effect of cross winds, atmospheric dispersions along the flightpath, and whether 
there is active guidance during any phase of the flight to correct accumulating errors. 

For some flight projects the vehicle may reach the target area at supersonic speed. In this case all data 
must be transmitted and recorded before impact since the vehicle may be destroyed on impact. 
Sounding rockets often rely on parachute systems for terminal descent and a beacon for the recovery 
of payload or experiment capsules. The parachute also makes it easier to track, locate, and recover 
the data packages. 

Parachute deceleration may be desirable for SRV capsule recovery. A parachute, even a small one, 
has the following advantages: easier tracking; more certainty of quickly finding the entry capsule; 
and a large decrease in g-loads at impact. The disadvantages include: increased system complexity 
associated with parachute deployment; a slight increase in the probe mass; and perhaps a lower static 
stability margin during hypersonic flight. Many of these questions and concerns can be resolved with 
sounding rocket flight tests. 
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8. DESIGN PROCESS SUMMARY 


NASA 's Stardust sample return mission returned safely to Earth when the capsule carrying cometary 
and interstellar particles successfully touched down at 2:10 a.m. Pacific time in the desert salt fiats 
of the U.S. Air Force Utah Test and Training Range. 

Ten years of planning and seven years of fight operations were realized early this morning when we 
successfully picked up our return capsule off of the desert floor in Utah. 

The Stardust spacecraft released its sample return capsule at 9:57 p.m. Pacific time last night. The 
capsule entered the atmosphere four hours later. The drogue and main parachutes deployed at 2:00 
and 2:05 a.m. Pacific time, respectively. The sample return capsule's science canister contained 
comet and interstellar dust. NASA 's Stardust spacecraft traveled 2.88 billion miles during its seven- 
year round-trip odyssey. 


NASA JPL Announcement, January 15, 2006 


8.1 Typical Sample Return Mission (Stardust) 

Stardust is a NASA comet rendezvous and sample return mission. The Stardust sample return capsule 
successfully returned to Earth in January 2006. The spacecraft captured particles from the tail of a 
comet. Two parachutes were employed during entry; a small one at 35.8 km altitude and the main 
chute at 3.05 km altitude. The Stardust sample return capsule is a blunt entry body (60-deg sphere- 
cone) with a PICA heatshield and SLA-561 ablation materials to protect the interior from entry 
heating. The sample return capsule entered the atmosphere at a speed of 12.5 km/s (relative to the 
atmosphere) and a flightpath angle of -8.4 deg. (The Stardust SRV is called a sample return capsule.) 
This was the fastest Earth entry ever for an artificial object. 

Stardust is the first U.S. mission designed to return samples from another body since the Apollo 
missions to the moon; it is the first U.S. mission designed to obtain samples robotically in deep space 
and return them to Earth. Data for the successful Stardust mission are presented below to show 
typical values that could be encountered in a Mars sample return mission because the necessary entry 
velocity may be comparable (Jet Propulsion Laboratory Stardust Mission, 2006): 

Launch Vehicle: Delta II (3 stages plus 4 solid rocket boosters attached to first stage) 

Launch Date: February 1999 from Cape Canaveral 

Distance of travel (entire mission): 4.63 billion km 

Duration of Mission: 83 months in space 

Power: Solar panels (6.6 m ) 

Earth entry velocity (inertial): 12.8 km/s 

Landing site: U.S. Air Force Utah Test and Training Range 
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Size: 

Spacecraft: 1.7 m high x 0.66 m wide x 0.66 m deep 
Sample return capsule: 0.8 m dia x 0.5 m high 


Mass: 

Spacecraft 254 kg 

Sample return capsule 46 kg 

Fuel 85 kg 

Total 385 kg 


Program Cost: $168.4 million (not including launch vehicle) 

The Stardust Earth-entry trajectory was simulated using the NASA TRAJ trajectory simulation code. 
Body geometry and initial conditions are given in table 19. Results of the trajectory calculation are 
shown in table 20. Trajectory variables (velocity, deceleration, heat flux and total heat load) are 
shown as functions of altitude and/or time from entry in figures 30 through 34. Table 20 shows that 
peak heat flux occurs prior to peak deceleration and peak dynamic pressure. Table 20 also shows that 
less than 10 percent of the total heat load is due to radiative heating. Figures 3 1 and 32 show the 
altitude of maximum deceleration and the altitude of maximum total (convective plus radiative) 
stagnation point heat flux, respectively. 

Note the differences between inertial and relative velocities and between inertial and relative 
flightpath angles in table 19. The inertial values are the variables used in the equations of motion. 

The relative values are used to determine aerodynamic forces and moments and heat fluxes. 

TABLE 19. STARDUST TRAJECTORY SIMULATION INPUT DATA 


Variable 

Value 

Vehicle mass, m 

41.5 kg 

Cone half-angle, 9c 

60 deg 

Nose radius, r N 

0.229 m 

Base radius, r B 

0.406 m 

Corner radius, r c 

0.019 m 

Surface area, Sw 

0.619 m 2 

Ballistic Parameter, BP a 

52.6 kg/m 2 

Conditions at Entrv Interface 


Altitude at entry, y E 

135 km 

Angle of attack at entry, a 

0 deg 

Inertial velocity at entry 

12.8 km/s 

Relative velocity at entry 

12.456 km/s 

Inertial entry angle 

-8.20 deg 

Relative entry angle 

-8.428 deg 


a For altitudes greater than 37.5 km 
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TABLE 20. STARDUST TRAJECTORY SIMULATION OUTPUT DATA 

(NASA TRAJ CODE, 1999) 


a. Maximum values and corresponding trajectory points 


Variable 

Maximum 

Time, 

Altitude, 

Velocity, 

Value 

s 

km 

km/s 

Convective Heat Flux 

605.4 W/cm 2 

49.8 

64.1 

10.9 

Radiant Heat Flux 

102.9 W/cm 2 

48.3 

65.5 

11.2 

Total Heat Flux 

705.1 W/cm 2 

49.0 

64.8 

11.1 

Stagnation Pressure 

33,300 Pa 

60.0 

56.0 

8.1 

Dynamic Pressure 

16,600 Pa 

60.0 

56.0 

8.1 

Deceleration Magnitude 

315.7 m/s 2 

60.0 

56.0 

8.1 

b. Heat load at stagnation point 

Heat Load 



Value 


Convective heat load, Qconv 

18,290 J/cm 2 


Radiative heat load, Qrad 

1270 J/cm 2 


Total heat load, Qtotal 


19,560 J/cm 2 



Notes: 

1976 U. S. Standard Atmosphere 
Gravitational model includes J 2 harmonic 
Destination planet is rotating Earth 


8.2 Design Process for Sample Return Missions 

The analytical, computational and experimental methods described in this document can be used in 
the design of an entry capsule for a Mars sample return mission. The process has several stages; some 
stages may require iteration. 

8.2.1 Mission Design Requirements 

The process can begin with maximum or minimum payload weight requirements. The Earth entry 
velocity is determined by the mission profile. For a return from Mars the entry velocity could be as 
high as 15 km/s. Mission requirements may limit the temperature in the capsule interior in order to 
protect the sample. The sample must survive landing impact. 
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8.2.2 Initial Probe Design 


An overall conceptual design is developed based on initial assumptions for the fully deployed 
spacecraft along with SRV mass and stability considerations. The ballistic parameter is estimated 
from size, weight and estimated drag coefficient. This value of ballistic parameter is used to calculate 
atmospheric trajectories and perfonn parametric analyses of SRV variables. 

The key variables in vehicle-trajectory optimization are entry velocity, entry angle and ballistic 
parameter for nonlifting, ballistic entry vehicles. The effect of varying the entry angle and ballistic 
parameter are discussed in section 4. 

It is often possible to reduce maximum convective heat flux (and TPS weight) by reducing ballistic 
parameter and increasing flightpath angle at entry (more shallow entry). However, the total heat load 
may increase (requiring a thicker TPS) due to longer flight time. 

Initial estimates of convective and radiative heat fluxes can be obtained from the analytical 
expressions given in Section 4 or from trajectory simulation codes. Note that estimates of radiative 
heat flux are less accurate than estimates of convective heat flux and total heat flux estimates are only 
approximate. With these heat flux estimates, an approximate TPS thickness can be detennined and an 
approximate vehicle mass can be calculated. 

What ultimately drives the selection of a TPS material and its thickness distribution is bondline, or 
attachment, temperature and total heat load. Consideration of total heat load is extremely critical 
when samples are very sensitive to temperature. 

At later stages, more detailed flow field calculations are made to improve the convective and 
radiative heat flux predictions. This will allow more accurate determination of the TPS thickness 
distribution. 

8.2.3 Terminal Descent and Recovery 

One of the most important considerations in entry vehicle design is whether a parachute system 
should be included (see sec. 5). The advantages of parachute recovery are: a) entry vehicle tracking at 
high altitude is much easier; b) the probe or canister has a much lower ground impact velocity; and 
c) recovery operations are much easier (the probe is easier to find). Disadvantages may include: 
a) slight increase in mass; b) increased mechanical complexity with the separation of the probe from 
the aft-cover and parachute; and c) slight decrease in static stability margin. 

A small parachute may be able to adequately decelerate the entry vehicle. Even though impact-area 
drift may increase significantly, ground radar should be able to track the probe and parachute and 
ensure efficient recovery operations. 
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8.3 Systems Engineering and Life Cycle Considerations 

Systems engineering is the art and science of developing an operable system capable of meeting 
mission requirements within imposed constraints including (but not restricted to) mass, cost and 
schedule. 


Griffin and French, 2004 

Sample return vehicle design is just one element of the spacecraft mission life cycle process. 

Certainly SRV design and TPS material selection are critical elements and several iterations may be 
required to satisfy mission requirements and cost and schedule constraints. It is important to note 
some of the many other factors that are required for a successful spacecraft mission. NASA requires 
well-defined technical workflow and management reviews. These items are combined in NASA's 
systems engineering and technical integration process. NASA and many other organizations follow 
similar processes developed over many decades. NASA requires all significant projects to travel 
along this path from initial concept to mission completion (Shishko, 1995). 

Project management and program control are part of the NASA project cycle. NASA programs 
follow a sequence that consists of distinct phases separated by “control gates.” Failure to pass a 
control-gate review can result in project tennination (see table 21). 


TABLE 21. PROJECT LIFE CYCLE PHASE DESCRIPTIONS AND CONTROL GATES 


Phase 

Description 

Control Gate 

“Zero” 

Advanced studies, mission feasibility 

Mission Concept Review 

A 

Preliminary analysis, mission definition, 

Mission Definition Review 

B 

System definition, preliminary design 

Nonadvocate Review 
Preliminary Design Review 

C 

Design 

Critical Design Review 

D 

Development, fabrication and integration, 
assembly, test 

Operational Readiness Review 
Flight Readiness Review 

E 

Operations 

— 


Final consideration must be given to the many systems, subsystems, and support systems that must 
be factored into the project plan and schedule; these include: 

Launch vehicle, environmental control 

Ground support, data processing 

Computers, hardware and software 

Payload, lander and SRV 

Attitude and thermal control 
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Electronics, power, batteries, propulsion 
Communications, uplink and downlink 
Vehicle and component structure 

Operations support, tracking, command and data handling 

The spacecraft must be maintained on optimal or near-optimal trajectories from launch to orbit, 
interplanetary, descent, landing, ascent, return, and Earth entry. All of this must be accomplished 
without compromising mission objectives. 


8.4 Space Environment 

There are many additional concerns in the design of a spacecraft mission (Griffin and French, 2004 
and Fortescue et ah, 2003). A few of these related to the space environment are listed below: 

Vacuum - The near vacuum conditions in space lead to significant outgassing of materials. This 
outgassing can in turn result in changes in the electrical and mechanical properties of spacecraft 
materials. 

Plasma effects - The effect of dissociation and ionization of gasses around the spacecraft and the 
SRV must be considered (Van Allen radiation belt effects, bombardment of the spacecraft by 
electrons, electromagnetic interference, and communications blackout during Earth entry). 

Earth’s magnetic field - The primary effect of variations in the Earth’s magnetic field is on the 
spacecraft attitude control system. 

Other radiation effects - Solar flares, x-rays, and ultra-violet radiation can affect electronic systems. 

Planetary environment - Wind blown dust on the Martian surface as well as temperature extremes 
can affect the lander vehicle and the sample collection process. (Temperature extremes are also 
encountered in the interplanetary trajectories.) 

Meteoroids and orbital debris - Meteoroids and orbital debris can cause damage to the spacecraft 
and the SRV. 

Microgravity and weightlessness - Microgravity and weightlessness can affect spacecraft 
components and the flow of fluids. 
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APPENDIX 


Acronyms and Abbreviations 


AEDC 

AFE 

AFRSI 

AGARD 

AHF 

AIAA 

ANSI 

AOTV 

ASTM 

AVCOAT 

CFD 

CM 

DOF 

DSMC 

FIAT 

HAVOC 

ICBM 

IEEE 

IHF 

JPL 

LCA 

LEO 

NACA 

NASA 

NASP 

NIST 

OMS 

PAET 


Arnold Engineering Development Center 

Aeroassist Flight Experiment 

Advanced flexible reusable surface insulation 

Advisory Group for Aerospace Research and Development 

Aerodynamic Heating Facility 

American Institute of Aeronautics and Astronautics 

American National Standards Institute 

Aero-Assist Orbital Transfer Vehicle 

American Society for Testing and Materials 

Trade name for ablation material developed by AVCO Corp. 

Computational fluid dynamics 

Command Module (Apollo) 

Degrees of freedom 

Direct simulation Monte Carlo 

Fully implicit ablation and thennal response program 

Hypersonic vehicle optimization code 

Intercontinental ballistic missile 

Institute of Electrical and Electronic Engineers 

Interaction Heating Facility 

Jet Propulsion Laboratory 

Lightweight ceramic ablator 

Low Earth orbit 

National Advisory Committee for Aeronautics 

National Aeronautics and Space Administration 

National Aero-Space Plane 

National Institute of Standards and Technology 

Orbital maneuvering system 

Planetary Atmosphere Entry Technology 
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PICA 

Phenolic impregnated carbon ablator 

POST 

Program to optimize simulated trajectories 

PTE 

Panel Test Facility 

RCC 

Reinforced carbon-carbon composite 

RCG 

Reaction cured glass 

SI 

International System of Units 

SLA 

Superlight ablator 

SRV 

Sample return vehicle 

STS 

Space Transportation System 

TABI 

Tailorable advanced blanket insulation 

TPS 

Thermal protection system 

TPSX 

Ames database of TPS material properties 

TRAJ 

Trajectory code 

TUFI 

Toughened uni-piece fibrous insulation 


Nomenclature 

a angle of attack 

P inverse scale height 

X viscous interaction parameter 

8 boundary layer thickness 

s radiation emissivity of material surface 

(j) bank angle (measured from vertical) 

y flightpath angle (measured from horizontal, positive up); or isentropic exponent (ratio of 

specific heats, c p /c v ) 

X mean free path, or Earth latitude 

p dynamic, or absolute viscosity 

9 angle between tangent to body surface and freestream direction; or angle between body axis 
and body surface 

p mass density 

a Stefan-Boltzmann constant, 5.670 400 x 10 8 W/(m 2 -K 4 ) 

\j / heading angle (measured from entry plane) 

v kinematic viscosity 

t fineness ratio 
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co frequency of oscillation 

A axial force 

BP ballistic parameter 

c speed of sound 

C Chapman-Rubesin constant 

CG Center of mass (or center of gravity) location 

CP center of pressure location 

Ca axial force coefficient 

Cd drag coefficient 

C D() zero-lift drag coefficient 

Cf skin friction drag coefficient 

Cl lift coefficient 

C m pitching moment coefficient 

C pitching moment curve slope 

a 

C N nonnal force coefficient 

C N nonnal force curve slope 

a 

C p pressure coefficient 

c p specific heat at constant pressure 
c v specific heat at constant volume 
D drag force; or diameter 

g acceleration due to gravity 

h enthalpy per unit mass; or geopotential altitude 

I mass moment of inertia 

J2 zonal hannonic coefficient in a series representing the Earth's gravity field (primarily the 

effects of Earth oblateness) 

K arbitrary constant 

ki dynamic stability parameter 

k2 static stability parameter 

k thennal conductivity 

Kn Knudsen number 

L lift force, or vehicle reference length 

Le Lewis number 

m mass 
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M Mach number, or pitching moment 

M molecular weight 

N nonnal force 

p pressure 

Pr Prandtl number 

q dynamic pressure; or heat flux (quantity of heat per unit area and per unit time) 

Q total heat load 

r radius, or distance from Earth center to vehicle center of mass 

R gas constant, or Earth radius 
R Universal gas constant 

Re Reynolds number 

Reo Reynolds number based on diameter D, p V D / p 

Re L Reynolds number based on length L, p V L / p 

Re x Reynolds number based on distance x, p V x / p 

S reference area, or surface area 

St Stanton number 

T temperature 

t time measured from entry 

V freestream or vehicle velocity magnitude, or modified viscous interaction parameter 

W weight 

x axial, or longitudinal coordinate 

y geometric altitude, or coordinate normal to body surface 


Subscripts 

0 sea level conditions 

1, 2 conditions upstream and downstream of shock wave, respectively 

AW adiabatic wall conditions 

B conditions at base of body 

C cone, or corner between cone surface and base 

CG center of gravity 

COND conductive heat transfer 

CONY convective heat transfer 
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CP 

D 


e 

E 

F 

FLIGHT 

INCOMP 

L 

LAM 

M 

MAX 

N 

RAD 

RE-RAD 

SPACE 

t 

T 


TOTAL 

TRANS 

TURB 

W 

WIND TUNNEL 

x 

oo 


center of pressure 
diameter 

boundary layer outer edge conditions 
conditions at entry interface altitude 

characteristic fluid temperature (usually mean or freestream) 

flight conditions 

incompressible flow 

reference length 

laminar flow 

molecular scale 

maximum 

conditions at nose of body 

radiative heat transfer (from shock layer gases to vehicle surface) 
re -radiative heat transfer (from vehicle surface to space) 
radiation to space or radiation si nk 
stagnation or conditions at stagnation point 

approximate location where boundary layer transition (from laminar flow to 
turbulent flow) occurs 

convective plus radiative heat transfer, CONV plus RAD 
boundary layer transition 
turbulent flow 

wall conditions or body surface area 
wind tunnel conditions and model scale 
distance along surface flat plate 
freestream conditions 


Units 


atm 

atmosphere 

bar 

1CT 5 Pa 

Btu 

British thennal unit 

deg 

degree 

ft 

foot 
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g gram 

h hour 

in inch 

J joule 

K Kelvin 

lb pound mass 

lbf pound force 

m meter 

mi mile 

min minute (time) 

mol mole 

N newton 

°C degree Celsius 

°F degree Farenheit 

°R degree Rankine 

Pa pascal 

psf pound per square foot 

psi pound per square inch 

rad radian 

s second 

slug unit of mass in foot-pound-second system 

W watt 


Constants and Conversion Factors 

1 deg = (tt:/ 1 80) rad = 0.01745329 rad (plane angle) 
n = 3.14159265358979 ... 
e = 2.7 18281 828 . . . (base of natural logarithms) 

1 Btu = 1055.056 J (International Table) = 1054.350 J (Thermochemical) 
1 J = 1 N-m 
1 W = 1 J/s 
1 Pa = 1 N/m 2 
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Length 

1 in = 2.54 cm 

1 ft = 0.3048 m 

1 mi = 5280 ft = 1609.344 m 

1 nautical mile = 1852 m = 1.15077945 mi 

Area 

1 in 2 = 6.4516 cm 2 
1 ft 2 = 0.092903 m 2 

Volume 

1 in 3 = 16.387064 cm 3 = 0.016387 liter 
1 ft 3 = 0.028317 m 3 = 28.317 liter 

Speed 

1 mile per hour = 0.44704 m/s 
1 ft/s = 0.3048 m/s = 0.0003048 km/s 

Force 

1 lbf= 4.448222 N 

Pressure 

1 bar = 100 kPa = 100,000 N/m 2 
1 psi = 1 lbf/in 2 = 6894.7572 Pa = 6894.7572 N/m 2 
1 psf = 1 lbf/ft 2 = 47.88026 Pa = 47.88026 N/m 2 
1 atm = 101,325 Pa = 760 mm Hg = 760 torr 
1 atm = 1.01325 bar = 14.696 lbf/in 2 

Mass and Density 

1 lb = 0.45359237 kg 

1 slug = 14.59390 kg 

1 lb/ft 3 = 16.018046 kg/m 3 = 0.031081 slug/ft 3 
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Energy and Power 

1 ft-lbf = 1.355818 J 

1 ft-lbf/s = 1.355818 W 
1 Btu/s = 1055.056 J/s = 1055.056 W 
1 Btu/(ft 2 -s) = 1.135653 x 10 4 W/m 2 (heat flux) 

1 Btu/°R= 1899.101 J/K(heat capacity) 

1 Btu/(lb-°F) = 4186.8 J/(kg-K) (specific heat capacity) 

1 Btu/lb = 2.3218 x 103 N-m/kg (specific enthalpy) 

Fluid Properties 

1 lb/ft-s = 1.488164 Pa-s (viscosity) 

1 Btu-in/(s-ft 2 -°F) = 519.2204 W/(m-K) (thennal conductivity) 

1 Btu/(s-ft 2 -°F) = 20,441.75 W/(m 2 -K) (heat transfer coefficient) 

Temperature 

T(K) = 5/9 T(°R) 

T(K) = T(°C) + 273.15 
T(K) = (5/9) [T(°F) + 459.67] 

T(°C) =(5/9) | T(°F) -.32 1 
T(°R) = T(°F) + 459.67 
1 °C = 5/9 °F 

International System of Units (SI) Prefixes 


n 

nano 

1 0 9 

fi 

micro 

10“ 6 

m 

milli 

10 3 

c 

centi 

10 2 

k 

kilo 

10 3 

M 

mega 

10 6 

G 

giga 

10 9 
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Figure 1. Characteristics of hypersonic flow around a blunt body. 
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Figure 2. Shock layer temperature for equilibrium air (Gordon and McBride, 1994). 



Figure 3. Vehicle flight profiles in Earth atmosphere (Hansen and Heims, 1958 and Howe, 1990). 
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Reynolds number, Re 

Figure 4. Flow regimes based on Mach number and Reynolds number. 
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Figure 5. Ratio of specific heats vs. temperature for air. 



a. Normal shock wave 


b. Oblique shock wave 


Figure 6. Normal and oblique shock waves. 
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Figure 7. Boundary layer on a flat plate. 
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Figure 8. Skin friction drag coefficient vs. Reynolds number for flat plate in incompressible flow. 
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Figure 9. Temperature vs. altitude, 1976 U.S. Standard Atmosphere. 
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Figure 10. Atmospheric density vs. altitude. 
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Figure 11. Mean free path vs. altitude, 1976 U.S. Standard Atmosphere. 
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a = angle of attack 
M = Pitching Moment 
M = - N x(CP-CG) 

CG = Center of Gravity 
CP = Center of Pressure 

Figure 12. Aerodynamic forces and pitching moment. 


A = Axial Force 
L = Lift 
D = Drag 

R=N+A=L+D 
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Figure 13. Hypersonic lift and drag for a flat plate at angle of attack. 
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Figure 14. Hypersonic drag coefficient for sphere-cones. 
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a. Sharp-nosed body with attached shock wave 



b. Sharp-nosed body with detached shock wave 


Figure 15. Sharp-nosed bodies with attached and detached shock waves. 
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Figure 19. Space Shuttle Orbiter re-entry trajectory. 
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Figure 20. Maximum deceleration for steep, nonlifting entry. 
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Figure 21. Altitude at maximum deceleration for steep, nonlifting entry. 
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Figure 22. Hypersonic pitching moment curve slope for sharp cones. 
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Figure 23. Peak stagnation point heat flux vs. entry velocity for steep Stardust entry. 
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Figure 24. Peak stagnation point total heat load vs. entry velocity for steep Stardust entry. 
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Figure 25. Schematic diagram of charring ablation process. 
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Figure 26. Atmospheric entry trajectory and ablation constraint. 
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Figure 27. Effect of emissivity on ablation constraint for an axially-symmetric body. 
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Figure 28. Effect of ablation temperature on ablation constraint for an axially-symmetric body. 
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Figure 29. Effect of nose radius on ablation constraint for an axially-symmetric body. 
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Figure 30. Stardust altitude vs. vehicle speed (NASA TRAJ Code, 1999). 
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Figure 31. Stardust altitude vs. deceleration (NASA TRAJ Code, 1999). 
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Figure 32. Stardust altitude vs. stagnation point total heat flux (NASA TRAJ Code, 1999). 
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Figure 33. Stardust stagnation point heat flux vs. time from entry (NASA TRAJ Code, 1999). 
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Figure 34. Stardust stagnation point total heat load vs. time from entry (NASA TRAJ Code, 1999). 
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